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АВ5ТААСТ 


The aim of this dissertation is to present a new method of 
engineering analysis and design for complex control systems, This 
method is the time domain infinite matrix method, The formulation of 
the infinite matrix follows from the convolution summation of sampled 
data systems, The mathematical basis of the time domain matrix for- 
mulation is presented in a discussion of the applicable concepts of 
infinite matrices and sequence spaces, This method of analysis and 
design is applicable to both continuous data and sampled data systens, 
For continuous systems it is necessary to introduce a fictitious sam- 
pler and hold of sufficient sampling rate to effect an accurate ap= 
proximation, 

The time domain matrix method is presented and illustrated as a 
method of analysis and design of linear, nonlinear, and time varying 
systems of the continuous or sampled data class, Sampled date, time 
varying systems may not be conveniently investigated by any other 
existing method, Furthermore the investigation of nonlinear systems 
is greatly simplified by the time domain approach, Multiloop sys- 
tems may be treated with ease and the signals at intermediate points 
throughout the loops are readily available, Also, systems with 
multiple nonlinearities may be investigated, for which there is not 
a presently available method of analysis and design, 

Two methods of design of a discrete compensator for a sampled 
data system are presented, These methods are accomplished directly 
in the time domain and allow for a compromise of specifications in 


the time domain, Also the response between sampling instants is 


71 





accounted for in one of the two design procedures, 
The time domain matrix method may be readily programmed on a 
digital computer and therefore provides a rapid analysis and design 


technique, 
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CHAPTER 1 
INTRODUCTION 


During the last two decades, feedback control systems have become 
increasingly important to our technolcgical civilization, and have par- 
ticularly contributed to national defense, Simultaneously, there has 
been increased interest and effort in the investigation of automatic 
control systems with respect to their analysis and design with various 
excitation signals, 

The I,R,E, defines a feedback control system as a control system 
comprising оспе or more feedback control loops, which combines functions 
of the controlled signals with functions of the commands to tend to 
maintain prescribed relationships between the commands and the con- 
trolled signals, Feedback control systems are a large class of sys- 
tems which include many subclasses of which a few are linear systems, 
nonlinear systems, multivariable systems, time varying systems, sampled- 
data systems, and adaptive systems, The active and dynamic elements 
are not limited and may be electronic, zlectromechanical, hydraulic, 
or pneumatic, 

An important subclass of control systems, sampled-data control 
systems, is a dynamical system which operates with sampled or quantized 
information, That is, the information is present as & sequence of dis- 
crete numbers in time, in contrast to continuous data systems for which 
the controlling information is monitored continuously in time,  Ordi- 
narily, the information is carried in the amplitude of the samples and 
may be considered pulse amplitude modulation, The block diagram of a 
typical sampled-data control system is shown in figure l-1, where the 
sampled-controller may be a special or general purpose digital computer, 
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Figure 1-1, Sampled Data Control System Block-Diagram 











The wide use of digital computers today has further stimulated inves- 
tigation of sampled-data systems, The seemingly unlimited possibilities 
of the use of a digital computer as an active element allow the possi- 
bility of considering these systems as a beginning towards the realiza= 
tion of the much discussed adaptive and learning control systems which 
may imitate the human brain and nervous system, 

There exist at present three major methods for the analysis and 
design of sampled-data systems as follows, 

l, Difference equation approach 

2, Frequency response methods 

3. 4 transform and modified transforms and the use of the root 

locus, 
All these methods assume that sampling occurs instantaneously or there- 
fore, that the pulses are of negligible width, The first method uses 
classical methods of difference equations and yields a solution only at 
the sampling instants, Frequency response methods are an attempt to 
extend concepts from continuous systems and suffer from the attendant 
disadvantages, The z transform uses a complex-variable transformation 
and determines the performance of the system by location of roots on 
loci in the complex z plane, and by inverse transformation, 

It is the purpose of this dissertation to present a method of en- 
gineering analysis and design for complex control systems, This method 
is the use of infinite matrices in the time domain, Analysis and de- 
sign for this method takes place directly in the time domain and avoids 
the necessity of transformations in complex-variables, This method has 
proven to be potentially useful for the investigation of a variety of 


systems, 





Time domain analysis and design has the important advantage of 
affording investigation directly in the domain of interest and direct 
evaluation of performance, therefore avoiding use of correlation 
theorems which are complex and may be inaccurate, This method avoids 
the difficulty of solution that is present for higher order systems 
using the z-transformation, In fact, the amount of work necessary for 
investigation is approximately the same for a first order as for any 
е order system, 

This method can be applied to many classes of systems which either 
do not lend themselves or are not possibly investigated by frequency or 
z-transform techniques, Therefore, beyond linear system investigation, 
one may analyze and design nonlinear and time-varying systems, Also 
the investigation of continuous systems is made possible through the 
introduction of a fictitious sampler and hold of suitably high sampling 
rate, The error introduced by this approximation of a continuous sys- 
tem by a sampled-data system can be reduced to a negligible amount with 
an attendant increase in labor of calculation, 

This dissertation is presented in three parts, The first part 
comprising Chapters 2 and 3 presents the mathematical background and 
formulation of the time domain matrix equations, The second part com- 
prised of Chapters 4 and 5 presents the analysis and design methods and 
verification of the theory by application to various types of systems, 
The third part comprised of Chapters 6, 7, and 8 is concerned with the 
design and realizibility of the digital compensator in the time domain, 


and the final conclusions and possibilities for future investigation, 





CHAPTER 2 
THE TIME DOMAIN MATRIX 

2-1 Introduction 

The method of investigation of sampled-data control systems most 
comnonly used today is the z-transformation, as is correspondingly the 
s-transform for continuous systems, The z-transformation converts the 
difference equations to a set of simultaneous algebraic equations which 
may be solved for the unknown variable and then by means of the inverse 
transform yields the response in the time domain, Formulation of the 
problem directly in the time domain allows one to avoid this transfor- 


mation and inverse transformation, 


2-2 A Sampled-data System and the z Transform 


A block diagram of a simple open loop sampled-data system is shown 





in figure 2-1, The definition of the z-transform of x (t) is E 
со 
E m i -n 

Ale) = Ran Ing) = x(nT)z (2-1) 

n=0 
EST _ ; 
where z =e = Uu + jv 
and T = sampling period 


Then it can be shown that the output at the sampling instants Y(z) is? 


Y(z) = G(z) X(z) (2-2) 
where G(z) is the z transform of G(s), 
In order to obtain the response between the sampling instants 
Jury“ introduced the modified z-transform where the response between 
the samples is; 


Y(z,n) 2 G(z,m) x(z) (2-3) 





X(s) 3 (s) G(s) Y(s) 


Syncronized Samplers T м 
with period T ЕС ---776) 
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Figure 2-1, Open Loop Sampied System 


Figure 2-2, Ап Open Loop Continuous System 


ON 


In order to obtain the response in the time domain, the inverse 


transformation must be utilized of which one formis, 


y(nT) = — $ ER au dz (2-4) 
Е 


where [' is the path of integration in the 2-р1апе 
which encloses all the singularities of the integrand, 
It is found that the usefulness of the transform lies with the use of 
the root loci on the z-plane, However, one desires to obtain the res- 
ponse directly in the time domain and that will be the subject of the 


next section, 


2-3 The Convolution Sumnation 
The time response of a continuous system G(s) as shown in figure 
2-2 to a continuous input x(t) is given by the convolution integral? 


t 


y(t) = қ glt-A) x(A) dA (2-5) 
= 00 


where A = dummy variable time delay 
апа g(t=A) = the delayed impulse response of the system G(s), 
This method of obtaining the response is usually avoided in favor of the 
Laplace transform due to the difficulty of evaluation of the integral, 
However, for the sampled data system as shown in figure 2-1, the 


input signal may be written; 


Xy 
Со 
ж X2 
x (t) = x(nT) 5(t-nT) - d 
2 (2-6) 
n=0 9 
where [x] is a column E 
matrix of n elements En 


where n extends to infinity 
y 


Now, in contrast to equation 2-5, one may write an equation for the con- 


tinuous output response as a convolution summation, 


п 


к= > gnT -kT) x(kT) (2-7) 
k=0 
n n 
or y pa 1 ж шо, 98 у = Š ВХ (2-8) 
k=0 k=0 


where Bn-k is often referred to as the weighting sequence, 
The values at the sampling instants of the impulse response, Es are 
related to g(t) the impulse response as; 
E ER = (2=9) 
t = nT 


Therefore, when the input is an ideal impulse of unity height, then; 


1 1 fork=n 


x 
ко 46 k#n 


and B LES the weighting sequence, 


If equation 2-8 is expanded, one obtains; 


п= 0 Yo = Ео-о Хо 
Bil Ya " So Xo * Go Xy 
a 2 2 У2 = 62-0 Xo + 8; ә * E2-2 X; 
о о о (2-10) 
п-п Yn 7 Eno Xo * Еп Xy + Ena XK + °? è o 


4 
£n-n *n 


And equation 2-10 may be written as? 
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It can then be 


Yo 


Уз 


matrix form as? 


Yo 

Уі 

J2 

J3 

Jn 
or 
where 
and 
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X | are column matrices of order r 


tE g ed 


m 8n-1 Xy T 
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(2~11) 


seen that this set of equations 2-11 may be written in 


(2-12) 


is a square matrix of order n called the system trans= 


fer u 


For the case of a time-varying system one then obtains instead of equa- 


tion 2-7 and 2-8? 


and the transfer matrix form for the system G(s) would be? 


y(nT) = > g(nT,kT) x(kT) 


n 


k=0 


(2-13) 


g(0,0) 0 0 0 — 0 


g(1,0) g(1,1) 0 0 0 
[с] 511770) E) (2729) 0 О 
; Я : 1 . (о-14) 


© ° ° © o 


o o 


gn,0) g(n,l) в(п,2) gln,3) » +» >» g(n3) 
For example,in order to obtain the output response for a unit 

step input to a system which is simply an integrator, we have to de- 

termine the X апа б matrices, For a unit step input, the input at 


the sampling instants is always one or, 


JE 


To determine the system matrix one must determine the weighting 


eH HH 


HF o 


values at the sampling instants, If the period T is equal to one second 
then one simple method of determining the sequence is to find G(z) and 
to divide; 





G(z) = 2 =] + fe + “i + 279 E e (2-15) 


Another, more generally useful method is to determine the impulse 
response g(t) and find the values at the sampling instants, For an іп- 
tegrator, G(s) = 1/s or: 

g(t) = u(t) a unit step, 


Therefore, the system matrix is. 


ENS 0 70 «0 O0 
TON 0 0 

GE EE So TTN 
л i Sona 


10 


= 


Then, the output can be obtained using matrix multiplication:” 


Y] = [e] x] 


1 0 1 1 
Ill 2 

d SENE > | (22107) 
7 


o 
o 9 
d 
0 


© 
0 


з 
H 
- 
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2-4 The Evaluation of the System Transfer Matrix of a Transfer Function 

In order to determine the output response of a system, the system 
transfer matrix must be available, This is determined with the most 
facility and accuracy by determining the g(t) by taking the inverse trans- 
form of G(s), In order to evaluate the matrix values substitute t = nT, 
This method of evaluation of the matrix is illustrated in Appendix A 
and Table А=1 gives the values for some representive systems, 

In many cases, a system has an undetermined G(s) or frequency res- 
ponse and the G(s) must be determined experimentally, In this case, it 
would be as convenient to determine the impulse response of the system 
directly and therefore the En values, For a large percentage of the 
sampled systems, a hold circuit filters the output of the sampler as 
shown in figure 2-3, In this case, the unit impulse is converted to a 
unit pulse of one sampling period width, Therefore, the z-transform 
equation for the output Y(z) is? 

Y (z) 2 GG(z)X(z) 
where G,G (2) 2 (5,6 (5 | (2-18) 
Therefore, we are interested in determining the matrix 
ЕЗ y and experímentally the values of this matrix may be 
determined by exciting the system G(s) with a unit pulse of period T, 


ШТ 





This method has been verified experimentally and yielded values within 
three percent of the expected elements of the matrix, This method of 
determination of the pulse response of the controlled system is easy to 
accomplish, and quite useful in investigating components with unknown 


transfer functions, 


2-5 Response Between the Sampling Instants 

In any but the most well behaved system, the response between the 
sampling instants is of interest and must be determined, For this pur= 
pose the modified z transform was introduced in the z domain, and an 
analogous method must be determined for the time domain matrix, Writing 
equation (2-9) one has; 


g, 5 Е(аТ) - g(t) (2-19) 
t = nT 


for the values at the sampling instants, Therefore, if the values at 
half-way between each sample are to be determined (m = 1/2 in the modi- 
fied transform) one has? 


g, (m = 1/72) = eat + ST) = elt) | 
ъ= пт + p (2-20) 


The symbol m was chosen to be consistent with the modified z-transform 
and is defined as the percentage of the period from the sample point, 
as shown in figure 2-4, The index m can assume a value O to 1 and is 
at the n sample when m - O and the желе sample when m = 1, 

For example, if the output between the sampling instants is required 
for the system shown in figure 2-5, then one writes the equation for the 
output as, 

Y (m) | = [6] x | (2-21) 


The impulse response of the system is; 
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х(%) ж 


G(s) 


Figure 2-3, Open Loop Sampled System With Hold 





(n-3) (n-2) (n-1) n 


Figure ?-4, The Intersample Response 
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g(t) = = 
Then, the values of the transfer matrix at the sampling instants 
(m=O) isn 

Ena В 


and the matrix is, 


1 0 0 
[ s l 0 ( | 
GI = 2-22) 
472 E 1 1 
-n 
e 


In order to determine the values at half-way between the sample 
points one must determine [6 ()] where m - 1/2, Therefore 
g, (1/2) = ea + 172) 


¿Y 2 0 0 m 
€ (в4/2)| m er (б (2-23) 
7245 „15 „-.5 


Then, the output at the sampling instants and at the mid-point of the 


sampling period is 


1 0 ON E 1 
J = | (.2679) i о 1 | = 1,3679 (2-24) 
Wa a) | 1,5032 
an 6065 0 0 Т „6065 
0/2] 2231 .6065 0 1 | = ,829 (2-25) 
‚0821 „2231 6065| 1 ,9117 





The determination of these values is discussed further in Appendix A 


and values are given for the systems considered, 


2-6 The Formulation of the Matrix Equation for Closed Loop Systems 

Before considering the closed loop system, one must consider the 
two block open loop system as shown in figure 2-6 and examine the mat- 
rix algebra, 


Then, the necessary equations are, 

B | = |с, | x| «м | = [о, |в] (2-26) 
Therefore, Y | = [ б, | | Gy | x| (2-27) 
and, in general, since matrix multiplication is not commutative it is 
incorrect to write the transfer matrices in the reverse order, that is, 


Y] = [ca] [c] x] ЕСІ [62] x] (2-28) 


Now, for a closed loop control system as shown in figure 2-7 one has; 
ж 4 
e(t) =r(t) - (c)t and e (t)sr (t)-e (t) 
In matrix form one obtains; 


Е| = R|- c] and c| = [5 [| E (2-29) 


Therefore, one has, 


E]= r|- [с | [9 | Е] ашпа E) = fr] + [e2][]} М R] (2-30) 
с] - [e] [e] {[21+ [ax] [old А] (2-31) 


Therefore 


1 о 87710 
@ 1 

where [т] = identity matrix = |O 0 so 0 
q 0 1 

апа [А] = inverse of matrix A, Thus, the solution for the output 


response inVolves matrix multiplication, addition, and inversion, 


Fortunately, one finds the inversion process is simplified by the fact 
"S 


х(%) 


Figure 2-5, 


Figure 2-6, 





First Order Sampled System 


Gs (s) 


Two Block Open Loop System 
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that the matrices are all lower trangular matrices; that is, all the 
elements above the diagonal are zero, Also, one finds equation 2-31 
may be written as? 

c | = [x] R| where [ K]= BB + [A J? (2-32) 
were — [a]- [az] (в| = the product of the transmission matrices 
in the forward path, 

Post multiplying equation 2-32 by [I] + [A] one obtains: 

27 А0) = [4 — 
Adding — [r]- [1] s [0] to both sides, then: 

[к]{[т]+[А]} # [а] +[:]- [:]#}][т]+[а]|-[:] 
Therefore, postmultiplying by the inverse of [I] + [A],one obtains: 

[к]= [т]-{[т]+ [А] (2-34) 
uhich substitutes subtraction for more difficult multiplication neces- 


sary in equation (2-32), 


2-7 The Solution for the Response of a Simple System 

In order to evaluate the response of a closed-loop system, the in= 
verse of the matrix [1)+[4]must be determined, Since a physical system 
always is represented by a lower triangular matrix, one powerful method 
for inversion is given by Frazer, Duncan, and Collar and presented іп 
Appendix в,“ Consider a simple approximate first order transfer func- 
tion system where the sampling period is one second and a) l, so that 
a neglible delay will be introduced, but there will be no output at tke 
п = 0 sampling instant, The system with a step input is shown in 
figure 2-8, 

The system matrix is determined in Appendix A and given in Table 
А-] and rewritten here? 
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О 0 0 о 9 о 0 
тоо 
ӘРЕ І 0.-% 0 
[2] =|1 1 : 
1 1 1 о 9 e 0 
then 
1 0 0 о е 
Ш +61 = 13 37 ага 
1 0 о 0 
_1 -1 1 0 0 
fir] “"ТОЛ/--| 0-1 1..0 
A 
Therefore, 0 0 1 
0 0 
т-тү [ = 10 
0 1 е о ө 0 
O | 
so that = — | 
0 O 0 1 0 | 
1 0 0 1 1 
E 01 0 |1|s5 1 
00 CES l | 
: : 1 | 
Кк» ӨТІ ғ 
Е. | 1. 


= 


(2-35) 


(2-36) 


(2-37) 


(2-38) 


The values between the sampling instants after the first period 


would be the same as at the sampling instants, since 


Gm = о) = (вв) Гог 121. 
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Figure 2-7, Closed Loop Sampled System 





Figure 2-8, Closed Loop System with A Step Input 
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СНАРТЕВ 3 
THE MATHEMATICAL THEORY OF INFINITE MATRICES 

3-1 Тһе Mathematical Theory 

It was shown in chapter two that the output time response could 
be obtained with the use of a matrix formulation, Furthermore, these 
matrices consist of elements whose values are those of a time response 
evaluated at discrete instants. For an open loop sampled-data system 
the matrix equation 2-12 will be rewritten here; 

y] = [c]x] (3-1) 

The column matrices x] and Х| contain the values of the discrete 
response and input respectively, If the response is to be determined 
for all time, then the order of the colum matrix, and the square sys- 
tem matrix, is infinite, ОҒ course, even if one does not need to evalu- 
ate the infinite number of response values, the system matrix can be 
considered as an infinite matrix, that is, of large order, The content 
of this chapter is a discussion of infinite matrices and infinite se- 
quences, and the calculation of the time response of a sampled-data 
system utilizing a matrix formulation in the time domain, 

The system matrix G is an infinite matrix since 

ГА| -А-(а,,) (1,1 1,2,3, 4-28...) (3-2) 

that is, the matrix A is an array of elements of infinite order,  Hence- 
forth, in this chapter, let [A] = A and the column matrices X | = x and 


5 


Y = y for convenience in notation, Cooke” discusses, at length, the 
characteristics of infinite matrices, Only the characteristics of 
immediate importance shall be mentioned here, In general, the theory 


of infinite matrices is connected with mathematical analysis and the 


theory of functions, 
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The operations of interest for the infinite matrices are, 


+ В = CT * ы. 


OO 
k=1 
M = y 


where à is a scalar, 

The input and output time responses are expressed as an infinite 
sequence of discrete values, These discrete values constitute the ele- 
ments of the column vectors x and y, A vector x may be considered a 
sequence space, A definition of sequence space is 1? 

A set S of sequences is called a sequence space when it contains 
the origin, and is such that, for every x and y in S and for every 
(complex) scalar c, x*y and cx are in S, 


The sequence space of interest is called g, the space of all se- 


quences, Then, the matrix equation 


у = Ах 
(3-4) 
n 
ше In = з? алк Ху 
k=0 


is a linear transformation of о on itself, that is, the system with 
matrix A transforms the input sequence into another sequence, the out- 
put sequence, When 


(3-5) 


“mk Ж anak 


the system is time invariant as previously discussed in section 2-3, 
The infinite, positive time, matrices for physical systems are al- 


ways lower triangular matrices (L,T,M,), that is? 
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813 = 0 when j? i, 
Also, for time invariant systems the elements are equal along the upper 
left to lower right diagonal, This matrix is called a diagonally in- 


variant matrix (D,I,M,) and exists when? 





as тар. for every і,1. (3-6) 


It is important to note, that a matrix multiplication is non-commu- 
tative for the class of systems which are time-varying (T.V.), These 
time-varying systems, for the operations defined in equation 3-3, are 
characterized by a non-Abelian (non-commutative) algebra, The time in- 

variant systems may be called a sub-class (T,I,) of the large class T.V, 
and are characterized by a Abelian algebra, That is, for time invariant 
systems, matrix multiplication is commutative, Thus for linear constant 
coefficient systems the order of the matrix multiplication may be re- 
versed, while for nonlinear or time varying systems they may not, 
Inversion of the infinite matrix A is often necessary and it is 
shown in reference 3, that if there is no i where ass = 0, and A is non- 
singular, then a unique inverse of A exists, The evaluation of the 


inverse of the L,T,M, is discussed in Appendix B, 


3-2 Convergence and Stability in the Sequence Space, 

In an automatic control system it is of great importance to deter- 
mine if the overall closed-loop system is stable and therefore the out= 
put time response is bounded, The output response is a sequence in the 
sequence space and it must be determined if the sequence is bounded, 
Therefore, if a step function test signal is applied to the system it 
is necessary to determine if the output sequence converges to a final 


value, If the system response diverges, the system is considered unstable, 
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For a closed loop control system, the output response may be written 
y = Ах (3-7) 
where A is the closed loop system matrix and 


A =0(1 +6)? (3-8) 


The G is an open loop system matrix as shown in figure 3-1 





Y(s) 





Figure 3-1. Closed Loop Sampled Data System 
Then, given an input signal of a convergent nature, that is, of a 
bounded nature, the problem is to determine if the output sequence y is 
of a bounded nature, Mathematically, the problem is, given a convergent 
sequence x, under the transformation A, does a convergent sequence y re- 
sult ? A theorem of fundamental importance concerning this problem is 
that of Kojima=Schur,? 
Kojima-Schur theorem; the necessary and sufficient condi- 


tions that 


у= Ах or у= 7  agy*X, (3-9) 
k=1 


should tend to a finite limit as n-» oo whenever X, is convergent are that 





со 

(а) > UT [Ê M for every n (3-10) 
к= 

(b) lim ай <ұ for every fixed k, (3-11) 
n> ب‎ В 
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GO 
(c) lim > a = <“ (3-12) 


nzeo k=] 
Moreover, if lim x, = x then the final value is 
k +00 
| 
co 
@) 7 = un y = eb +> «(x-Y (0-13) | 
n=» & =1 


A matrix satisfying a,b,c is called a k-matrix with e, and « for its 
characteristic numbers, For, the special condition of x = O in condi- 
tion (b), one defines the matrix A as a T matrix, These are the condi- 
tions for convergence and clearly define the stability of a control sys- 
tem in the time domain, By examination of the impulse response of stable 
systems a physical understanding of these conditions results, Consider, 


for an example, a time invariant open loop system with a transfer function 


G(s) = + and T = 1 as discussed in section 2-3, 


Then a, =a and a = 1,00 where m=ne-=k, It is obvious that cone 
nk n-k m 


dition (b) yields е = 1,00, Furthermore, it can be seen that this 


k 
system does not satisfy conditions (a) and (c). That is, examining condi- 


tion (a) one finds for n í 
oo 


that > | ak | = © e This result is as expected since the response 
K=1 


is diverging as was determined in section 2-3, equation 2-17, Therefore 


condition (a) requires that the series > [акі converge in order for 


the system to be stable, For a closed loop system of figure 3-1, it is 
| =] 
important at this point to recall that A= G (1 + G} 2 Therefore, А 
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тау satisfy the condition (a) while G may not, 
Theorem? If A is the matrix for a time-varying system (T,V,), the sys- 


tem is stable if and only if it isa Kr matrix, 


Corollary? If A is the matrix for a time invariant system, then it is 
a D.I.M. and A = tinj’ and A is stable if and only if Aisa 


T matrix, 
с 


For the large class of time invariant systems, the corollary implies 


that for stability it is necessary for 
Co 


m la, | < wheren=i-j (3-12) 


n=0 


For stable T.I. systems, one finds that «, = O, that is the matrix 


k 
is aT matrix, For condition (с) one obtains « = 1,00 for a type I 
servo system, Then, by the use of relation (à) one obtains, as expected, 
the final value as, 
SS = mK =X (3-15) 
The transformation on 0, accomplished by A as a ПЕ matrix is called 
a regular transformation, The property of importance is that as in 
equation 3-15, the T, matrices have the property of consistency, That 
is, every convergent sequence is transformec by such a matrix into another 
convergent sequence with the same limit when n approaches infinity, In 
addition, T, Or k matrices will frequently transform divergent sequences 
into convergent sequences, On this point it may be stated °° 
Corresponding to each unbounded divergent sequence i - and each 
bounded sequence Í у, { ; there is a general (square) T, matrix which 
carries ER. into 172 | 


In this chapter, the mathematical basis for the infinite matrices 
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has been discussed, and the operations defined, Furthermore, the theory 
of stability and convergence in the sequence space has been discussed, 
It is not to detract from this, that usually the conditions a), 5); с); 
d) are investigated simultaneously with the evaluation of the actual 


output sequence for a given systen, 








СНАРТЕК 5 
ANALYSIS OF CLOSED LOOP CONTROL SYSTEMS 
4-1 А Method of Evaluating the Response of a Closed Loop System Witn- 
out Inversion of Matrices 


Consider the simple, error sampled, closed loop system as discussed 


in chapter two and shown in figure 4-1, 


e (t) 





Figure 4-1, Error Sampled System 


It was shown that the output response time sequence may be written as, 


S -[el[r] «[e])^ 2] = [0-41 «lo? 2] 4 


Therefore, the evaluation of the output response involves the inversion 


of [I] + [G] as discussed in Appendix B. 

However, there is a simple method of evaluating the output response 
which avoids the inversion of a matrix, The output response column ma= 
trix may be written as, 

dq = [¢] £] (4-2) 
мһеге E] = the error sequence in time, 
The error sequence may be written as 

E] = Е] ~ С] (4-3) 
It can be seen that the error matrix can be easily evaluated by subtrac- 


tion, Then the error matrix multiplied by the system matrix will yield 
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the output response, Of course, the error matrix cannot be evaluated 
until the output sample values are available, However, it may be seen 
that the output sequence may be evaluated in a step by step procedure 
by examining the expanded form of equations 4-2 and 4-3 for the first 


three sampling instants, For a linear, time invariant system, one has, 


со Lo 0 O | eg 
cal =l 8 вк Ols (4-4) 
с; 82 81 &oļ| €2 

апа ео no са (ro = co) 
el = mnl- „wie (m - ce) (4-5) 
е; n Ca (г = c3) | 


Expanding equation 4-4 with row by row matrix multiplication one obtains 
the following results, 


For the first sampling instant, the time origin, one has, 


Co = 80 &o 


and (4-6) 


99 7 Yo 7 Co 
Therefore, solving for eg one obtains; 
ae = а. 


However, for all physical systems, there can not be an output immediately; 


that is, there is always a time delay of small, but real magnitude, There= 


fore, бо = О Гог systems of interest, Then one has 


wi 


ео = To 
and (4-8) 
со = 0 
For the second sampling instant one has? 


Ca 5 8460 % 8064 7 8460 = 81 Yo (4-9) 
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апа ы 
е, >= 0004 


Since g, and rg are known, then cy may be evaluated and then e, nay be 


evaluated, Now, for the third sampling instant one obtains; 


ti 


с; = 8260 % 8464 t 806; ^7 E2390 + 51е; 


апа (4-10) 
e5 7 Ta = Ca 


Therefore, it can be seen that one may proceed to evaluate the 
time response sample by sample in time, The equations such as 4-10 be= 


come lengthy for larger n and it is simpler to use the matrix form, The 


matrix form may be written, with gg = O, as, 


со 0 0 0 0 Го 

C4 = 84 0 0 (т, ез C4 ) (4-11 ) 
с 2 82 E4 0 © о о ne = б 2 ) 

^n En 8n-1 *n-2 (г, 7 е 


For a linear time invariant system, the system matrix is diagonally in= 
variant, Therefore, multiplication of | a] e] as a row-colunn multipli- 
cation may be replaced by a colum-colum multiplication, This is possi- 
ble since, for example, the third row is identical to the first column 


from the third element upwards, That is, one may write the a muitipii= 


cation as; 


eo ES 
cC = [Ж Eq 9 £o | e4 = Е бі» o е, (4-12) 
е; е; 
or alternately as; 
0 бе 
Co =| 81 eq |= 8290 + 844 (4-13) 
Ep 
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Multiplication is accomplished by multiplying the first element of 


E] by the last element of[G] of interest, 


Then, in general it is possible to write, 


0 ео 

| ба | ёч 
[c] = | 82| *2 
63| 9*3 

84| 94 


(45375 


If a table of the form of equation 4-14 is established, a step by step 


evaluation procedure may be used as follows; 


1) Evaluate eg = ro, and со = 0 


2) Then for c,, multiply the G and E matrices by starting at g, 


in the G mtrix, 


Then one obtains c} = g,€9. Now evaluate 


е, 2 Т) - с, апа fill inc, and e, in the computation table, 


at gz. Then one obtains c} = geo t £494. 


eg Zr, = C, and place the values in the table, 


3) Now evaluate c, by multiplying the G and E matrices by starting 


Then evaluate 


4) Continue this procedure for each sample point of interest, 


typically until the system settles to a final value, 


An example will best illustrate this procedure, 


matrix of figure 4-1 be; 


OOO 


е 


~OO00 


O O OTO 


Let the system 


(4=15) 


Assume a step input signal R | as о о o o 11 ә Then, es- 


tablishing a table as in equation 14 one calculates each value of e, 
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and e, step by step as outlined, This method is illustrated for the 
first two calculations by use of dashed lines to indicate the flow of 


the calculations as follows, 


G E C 
0 „= е = 1,0 бе = 0 
S > > 
2 pp 7 УСС = A S 
c] v 4ه‎ е; = Т, = C4 = d T. 40 = G0 — C4 = 240 (4-16) 
س س چ س ا‎ | 
„8 е; Ca 
1.0 еҙ Сз 


Then, the next step is illustrated schematically as, 


G E C 
o | Е БЕ н 
& 
С] = | деже = .60 с, = 40 | (4-17) 
77 ER 
86 е, =1-с, = - 040 | = | е, = 1.040 
ee Sl B 
1.0 | ез | | 8 


If this procedure is carried on step by step, one obtains for the 


response at the first six sample points, 


G E C 
0 1 0 
„% 6 wh 
EISE Ue | - хое (4-18) 
1,0 = 464 1.464 
1,0 ~ „3824, 1.3824 
LL 10 ~ „0358 1,0358 4 


To evaluate the intersample response one proceeds as-discussed in 


section 2-5, Rewriting equation 2-21 one has; 


са) | =| 0(m)] £] (4-19) 
If the G(m) matrix for the midpoint of the sampling period is found 


to bee 
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28 Y 0 О 0 
60 „20 0 О 0 
EJ = 90 ¿60 „20 0 0 (4-20) 
| 1,00 „90 ,60 „20 0 
Then one may calculate the intersampling response as, 
G (n E C 
20] 1 „200 
060 2600 2720 
ca/2) = | ,90 | -,040 |* 1.252 (4-21) 
1,00 = A464 1,423 
1,00 = „3825 1.209 


It is valuable to note at this point, that the calculations in the 
case considered in this section were simplified by the unity feedback 
condition which yields Е| = R] = cl. Single loop systems with other 


than unity feedback shall be considered in the next section, 


4-2 Evaluation b the Response of A Closed Loop System with Other than 
Unity Feedback, 
In this section the closed-loop system shall be considered with 
other than unity feedback as shown for one case in figure 4-2, The 
samplers are synchronized and with the same sampling period, Then the 


matrix equations may be written as? 


c| = [в | s| (4-22) 
P a al [ao] 
ТһегеҒоге Е | = R] = н| c] (4-23) 


Now, to evaluate the output response one uses equations 4-22 and 
4-23 and evaluates the Ca and Es step ty step as outlined in the previous 
section, The only difference is in this case, in order to evaluate El, 
one must evaluate [5] c | and then substract ba? tha sample value, from 


B to obtain the errcr e, at the particular sample, It must be pointed 
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Figure 4-2, Two Sampler Feedback System 


C (s) 





Figure 4-3, Error Sampled System 
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out that the solution in this case depended on the fact that the trans- 
fer blocks G(s) and H(s) vere separated by samplers, 
Consider the system shown in figure 4-3 which does not have samplers 


sevarating both transfer blocks, Then one may write; 


с| -|с|Е| (4-24) 
and E] = R]- В] (4-25) 
But, noting that B(s) = H(s)C(s) = H(s)G(s) Сй then; 
E| = r| -cl (4-26) 
z R|- [m] z] (4-27) 


Therefore, it can be seen that Е| is not readily available in equation 
4=26 since нс! is not usually available, 


Then, rearranging equation 4-27, one obtains ¢ 


| - » [н] | el = n] (4-28) 
-4 И [ва ee] (4-29) 


Therefore, in order to >. ا‎ the output response for this system, 


inversion of a matrix cannot be avoided, 


Finally, if one was analyzing a system such as shown in figure 4-4, 


the equations of interest are; 
El - R] - [nic] (4-30) 
ad Cc) - [a] [ol x| (4-31) 


If the response at an intermediate point in the system is desired, such 


as the output of G, (s), it is readily available by writing the following 


Ea] = 8| (4-32) 


The response at intermediate points in a system is of great importance 


equation; 


in nonlinear systems which will be discussed in section 4-5, 


34, 








Section 4-3 Analysis of Multiloop Control Systens, 

The Introduction of more than one feedback loop is necessary or 
inherently present in many control systems, The investigation of a 
multiloop sampled-data control system is complicated by the presence 
of samplers in some loops, and the absence of samplers in others, There- 
fore, all the feedback signals are not of the same form throughout the 
system, Consider at first a two loop system as shown in figure 4-5, 
which has samplers separating all the transfer blocks, 

One may write a set of simultaneous sampled-data equations, where 
the starred notation indicates a sampled signal 2 transform, as follows, 

E (s) -R (s) - E (sc (s) - M (s) 
M'(s) = EN (s)E (s) (4-33) 
C" (s) = бү (в)М (в) = G, (s)G. (s)E (5) 


Solving these equations simultaneously, one obtains for the sampled- 


output. 
a. (s)G" (s)g (s) 
с" (=) = AS E (4-34) 


l + G, (2) + б, (6/5: (в)н” (в) 


Тһе z-transformed equation then follows as? 


G, (z)G_ (z)R(z) 
E (4-35) 





Ci) = 
1 + G, (2) + G, (z)G, (2)H (2) 


One may then write the matrix equations directly from equations 4-34 
and 4-35, or alternatively one may derive the matrix equations directly 


from the system signals, In either case, one obtains? 


c| = A | 1] « [е |+ [c] [e | (ay \-1 e] | | (4-37) 
and C(m)]| = | 11+ [| +1 1111 үч” [а (а) |6,(а)| (4-38) 
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C (s) 





Figure 4-4, Three Sampler Single Loop System 





Figure 4-5, Multiloop Sampled System 
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If it was desired to avoid the inversion necessary in equations 4-37 and 
4-38, one may write a set of simultaneous equations which provide a means 
of step by step evaluation of intermediate signals, In this case one 


would write; 


E] = в]-м]- [Ric] (4-39) 
M|= [6] E] (4-40) 
с] = [|] м] (4-41) 


In order to evaluate El. M], and C | a step by step procedure is 
followed similar to that of the single loop method, First, assuming 
no immediate o'.tput for either Ga ог Gy ә one may write? 

my =O 
со 50 (4-42) 
and еу = го 
Then, for the next sample one obtains, 
IM = Big ео = Biga Yo 
where £,, = the first impulse response value of the Ga block and 


th 


= the n ” impulse response value of the Ga block, 


Ena 
Then, continuing, one obtains > 
Ca = gibo * Eop My ^ Ep =O (4-43) 
since рор = О апа т = 0, 
Therefore, one uses equations 4-39, 4-40, 4-41 in that order, step by 


step, For the next sample one obtains; 


eq = ту = mM = (hoe. + hco) = Ty - Era Yo - 0 
(4-44) 
since Cg = cg = 0 
and Шш» = Вга ёо * Бза ен (4-45) 
and finally 
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Сә 7 Bab Шо + Sab My = Esp My since mg = 0 

As is the usual case, it is actually easier and more methodical to use 
the matrix form throughout the calculations, Therefore, one will esta=- 
blish a table similar to that of section 4-1 using equations 4-39, 4-40, 
4-41, The flow of calculation is from E] to M] to с| and then back 
through the feedback loops again to E| , It is obvious, that actually 
one is simply following the signal around the closed loops, 

Finally, consider a multiloop system which possesses only one samp- 


ler in the control loops, The simple sampler typically is placed in the 


error channel as is shown in figure 4-6, 


C (s) 





Figure 4-6, А Muitilcop Sampled System 


One may write the equations for the signals as follows? 


х E. 
H- ENER 
М = GE 

A (4-46) 
C = GM =GGE 
В = HC 


Therefore, one obtains for the closed loop sampled error and output; 


Е = — T A (4-47 ) 








es? 


un ж 


апа X Ge R 





C= سفق‎ (4-42) 
1+ Go + HG G 


One may then write the matrix equations directly as, 
Е| ([ г] + [GÎ + [ос [} ч R] (4-49) 
с] GG, | f [1] +[0,] + (в) R] (4-50) 
The equation for the intersample response is then, 
C(m)| = [G,G, (m)] | Ме [д] +» ETA | аз) 


Therefore, whenever one desires to calculate the intersample response 


Bb 


and the inversion of the matrix I + Ga + HG G, has been already accom- 
plished, one simply evaluates the GG, (m) matrix and carríes out the 
multiplication, Therefore, multiloop systems with one sampler or samp- 
lers separating all transfer blocks may be equally treated by the use 

of time domain infinite matrices, The response at any sampler location 
is readily available and is usually of interest in the ínvestigation of 
multiloop systems, particularly with nonlinearities present, The intro- 


duction of nonlinearities into a control loop is treated in section 4-5, 


Section 4-4, Analysis of Time-varying Control Systems 

The analysis of time-varying sampled-data control systems may be 
accomplished using the time domain matrix method, A control system 
may have as one of the Өле function blocks, a component whose 
parameters are changing with time, This effect is present with high 
altitude jet aircraft, where the dynamic characteristics of the air- 
craft change with altitude, and therefore time, 

Consider the simple open loop system shown in figure 4-7, The 


impulse response of the time varying component is changing with time, 
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AA Cy Time-Varying c (t) 


Component G 
Figure 4-7, Sampled Open Loop System 


Therefore the matrix equation for the output may be written as in section 


2-3, equation 2-13 and rewritten here. 


б] = [stu] El (4-52) 
тз g(0,0) О О 
g(1,0)  g(1,11) 0 "a 


| 6 (m,x)| = g (2,0) g(2,1) g(2,2) ж SZ 


9 о 
9 9 e 


о e e | 
In order to substitute numbers in the G(n,k) matrix, the time variation 
of the system G must be known or determined, As an illustration, consi- 


der an abrupt change of 


G(s) = so from a = 1,0 to a = 2,0 at the third 


sampling instant (k = 2) where T = 1 second, Then, using the values of 


е. from table A-1 one has for the system matrix, 


0 0 0 0 0 04%: 

se) 0 О О О О 

7675 „3679 0 0 0 

„145 „7675 „вв 0 0 0 
[G(n,k)| = |.9685 us „тт 2898-00 0%": 


9884 9665. ,4960 44707 62838 шай 


о о 9 
o e e 
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r(t) 





Then, of course if this component had feedback introduced as shown 


in figure 4-@, one would have the following closed loop equation; 


ж 


с] = [G(n,x)| y En] + Сы R | (4-53) 


— >= = کے 





MR, Time Varying с(%) 
Component G 





Figure 4-8, Time Varying Closed Loop System 

In order to obtain the intersample response, the equation for the output 
response may be written, 

cm] =[c, „m| [ + [сок 2 $] (4-54) 
Analysis of a feedback system with more than one time varying element 
follows the same approach as for the time invariant systems, Te matrix 
equations are found to be the same as for time invariant systems with 
the time variation of a transfer function only affecting the system ma= 
trix itself, Therefore, for the system shown in figure 4-9, the follow- 


ing equation is obtained; 


a s : A 
с] = [e(n,k)| [н + [G(n,k)| lui] | R] (4-55) 
It can be seen that it is not necessary for samplers to separate the 
time varying component from all other transfer blocks, Therefore, if 


H (3) was time invariant, it would not be necessary to have a sampler 
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C (s) 
Time Varying 


G(s) 





ТЕ ee ER ee ds 


Time Varying 
H(s) 





Figure 4-9, System with Two Time Varying Components 
between G(s) and H(s) for the use of the time domain matrices, For 


the system shown in figure 4-10, one obtains the equation; 


С]: = [e(n] f EF] + [ан(а,к)) | - R] (4-56) 







Time-Varying 2 SR 
G(s) 


Figure 4-10, Time Varying System with A Feedback Component 


A control system, where one of the components is a time-varying 
amplifier is worthy of consideration, Consider the system as shown in 
figure 4-11, 

The gain a(t? is changing with time and therefore has a different 


magnitude at each sampling instant, The equation for the output may 
be written; 
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R (s) Time-Varying M(s ) C(s) 
A Ыр. Amplifier 4 G(s) 
A = a(t) 
Figure 4-11, Open Loop Time Varying System 
с| -(|мбыю) қ (4-57) 
ог 
Co 0 0 0 0 0 Го 
C4 2081 0 0 0 Do o I4 
C = | àgE2 a4 Es 0 0 0 гә (4-58) 
С» ад, 3848: 482834 0 0 г» 
5, ао0Е, %83 4282 4381 AS 
Then, the output at the third sampling instant is as expected, 
Co = a98aTo + а.г; (4-59) 
It can be seen that equation 4-58 may be written as follows, 
А 0 0 0 о а, 0 0 0 
c] =[a] la] r] =[g, o o o .|jo a 0 0 к] (4-60) 
Zo g& O 0 0 0 az0 
| 
l83 82 & O |10 0 O ag! 


the form of the matrices, 


sampler, 





tative as was previously discussed, 


The operation of multiplication of the A and G matrices is not commu- 


That this is so is obvious from 


Now, consider an open system identical to 
that of figure 4-11 except that a hold circuit immediately follows the 


Then, if the amplifier changes gain during the sample period; 
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this must be accounted for in the system matrix G, which is a function 

of m, In most cases it is reasonable to make the simplifying assump= 
tion that the gain changes only at the sampling instants, In other words, 
the assumption is that if the gain is changing continuously with tim, 

the time constant of this change is much greater than the sampling period, 
Therefore, one is approximating the gain by a staircase function as 


shown in figure 4-12, 





Figure 4-12, A Time Varying Gain 


Then, on this basis, one may analyze closed loop systems with a time 
varying gain, For the system shown in figure 4-13, one may write the 


equations for the system signals as; 


HI mM" [a] ce] (4-61) 


and C]= [c][a] E) (4-62) 
Equations 4-61 and 4-62 may be solved by the step by step procedure 


li 


i 


previously outlined for the time-invariant systems, Using this step by 
step method, the gain change at each interval is clearly displayed to 
the investigator, Alternately, one may solve for the output response 
and error response by means of inversion of the matrix in the following 


equations; 
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P(s) c(s) 


Time-Varying G(s) 


Gain A 


Figure 4-13, System with A Time-Varying Gain 


=] = (0211-8768) (40 > в] (4-63) 
at o] = [ol ГАС] [4]] =] (4-64) 


If there was no sampler between the output and the feedback block 


H(s), then one obtains the equation for the error sequence in the same 
manner as carried out in section 4-2, Therefore, one obtains for the 
error sequence; 
E] =(01] + Cao) (Ip? 8] (4-65) 

As an example, consider a system where H(s) = 1 as shown in figure 
4-13. Let R(s) = 1/s, and the gain function be a(t) = T Р, with a 
sampling period of one second, In this the time constant of the gain 
change is only double the sampling period, but the results are instruc- 
tive with this marked gain change as shown in figure 4-14, 


Use equations 4-61 and 4-62 where | H] = [т] апа one may write, 


E] = в] -0] (4-66) 
c] =[6][a]E] = [0] P| (4-67) 
vhere P] =[a] E] the amplified magnitude of the error pulse 


and 
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1,9 ешш 
Gain 
a(t) | ш 
s | 
.2- | l = | a 47, 3679 
| ! 78342,2231 
| | i i 
О | 
0 Т 2T 3T АТ time 
Figure 4-14. Time Varying Gain Example 
_ o Wo o 6 
cl Ш O 0 
Т oo 4. 0 
1,0 e oh MS 
LA? + j 
Then, one obtains; 
G A 
О 0 0 0 159 0 0 0 1 
۸ о о о О „6065 0 О 600 
Шз а м о о On Of ser „0544 
1,0 ЕЕ & 0 О О 0 22231 | -.2991 
boe. pU 
G P 
0 1.0 | 0 
А „ 3639 oh 
=| 8 „0200 | = „95% (4-68) 
1,0 == 06673 1,2991 
1.0 я 159552 


o 


The availability of the amplified magnitude of the error pulse as Р] 


is often useful іп the investigation of nonlinear systems which shall be 


discussed in the next section, 
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4-5, Analysis of Control Systems with Nonlinear Components, 

The analysis of sampled-data control systems with nonlinear com- 
ponents may be accomplished by the use of time-domain infinite matrices, 
The methods most commonly used for continuous nonlinear systems are the 
describing function and the phase-plane, These methods are somewhat 
limited in their application to sampled-data control systems, Consider 


a single loop nonlinear system as shown in figure 4-15, 






G(s) | € (8) 





Nonlinear 
Component 


R (s) e» | E (s) 





Figure 4-15, Single Loop Nonlinear System 
Then, it can be Seda that the nonlinear component has replaced the varia- 
ble gain amplifier of figure 4-13, The nonlinearities considered are 
those sensitive to the magnitude of the signal input; that is, they 
have a nonlinear amplitude response such as shown in figure 4-16 for a 


saturating amplifier, Е 







(Voltage ) 





in 


Figure 4-16, Saturating Amplifier Characteristic 
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Therefore, as the input voltage varies in magnitude, the gain of the 
amplifier changes, Since the input to the nonlinear element is the 

held error magnitude, this magnitude determines the gain, However, the 
error is evaluated directly in the time domain and therefore a gain cal- 
culated for a calculated error input is also known as the gain at the 
specific sample instant, Therefore, the nonlinear element may be * 
treated as a time-varying gain, This statement is true for hysteresis, 
a relay servo, and other single nonlinearities following a sampler, 

Then the output of the nonlinearity is the sequence P| o. One may write 


for the output response and error response, 


El = c] (4-69) 


c] =[a] [u] £] (4-70) 
where 
O 0 70.4. 
0 ш 0 0 
u] =| 0 O u Oo. <>. 
0 о 0 13 


The values of u, depend upon the error magnitude and are determined in 
the step by step solution, As an example, consider a saturating ampli- 
fier with a characteristic curve as shown in figure 4-17, which is the 
nonlinear component of figure 4-15. If the input signal is a unit step, 
then the amplifier is expected to saturate, This amplifier has a maxi- 
mum gain of two throughout the linear region, Then, using equation 
4-69, at the first sampling instant (n = 0) eg = 1,0 and therefore 


Up = 1.0, and pp = 1,0, a gain of one, The G matrix for this example is, 
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Figure 4-17, Saturating Amplifier Characteristic Example 


0 0 

А 0 

[с] = .8 4 

e „8 

|1,9 1,0 
Therefore, c4 7 g4eg = „40 and e, = .60, This magnitude continues 

to saturate the amplifier and p, = 1.0, and u, = 1.667, Continuing in 


this manner, one obtains the following matrix solution, 


c] = [a] [u] E] =[6) P)] A-A) 
0.. 1 0 0 0 0 0 1 
(ы... о 1,667 0 0 оо „© 
8 0 0 а ос = 
1,0 0 0 0 1,56 0 0 = 64 
1,0 0 0 Oo (Oo. %9 0 -.28 
1,0 0 0 0270202 + 424 
| 0 1 0 
4 1 ab 
S. 1.20 
Imo cal = iver 
1,0 | -.56 1,28 
1,0 | +.828 „576 
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The analysis of systems with a nonlinear component which is also 
a storage device follows along in the same manner, For example, if the 
nonlinearity was hysteresis in a magnetic amplifier, then the gain is 
dependent upon the past value of input voltage to determine which side 
of the hysteresis loop is applicable to the present signal, This can 
be seen clearly from figure 4-18 which shows a simple hysteresis loop, 
If the input was E, , there are two possible output magnitudes depending 
upon the magnitude of the input of the previous sample, If the previous 


sample had a magnitude of Eg, then the output magnitude will be Eye 





Figure 4-18, Hysteresis Characteristic 


Furthermore, any nonlinearity may be treated in this manner since the 
relation | 
Р] = ®] (4-72) 
applies to any nonlinearities whose characteristics are known or can be 
approximated, 
The analysis of systems with a nonlinear component which does not 
have a sampler immediately preceding it cannot be treated by this method, 
However, as an approximation a fictitious sampler may be placed before 


the nonlinearity, The accuracy of this approximation depends upon the 
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sampling rate and the time constants of the system, This approximation 
shall be discussed further in the next section, 

If there is more than one nonlinearity present in the system and 
samplers appear before every nonlinearity, then the investigation of the 
system by time domain matrices is entirely possible, In fact, for a 
system with many nonlinearities, analysis is possible with the introduc- 
tion of samplers with a high sampling rate, wherever a nonlinearity 
exists, The necessary sampling rate shall be discussed in the next sec- 
tion, 

Furthermore, the signal sequence magnitudes are available at inter- 
mediate points throughout the system and can aid an investigator in the 


analysis of the component requirements, Consider for example the system 


of figure 4-19 which has two nonlinearities, Then the equations for the 





Figure 4-19, System with Two Nonlinear Components 


error, the output sequence of G,, and the output sequence of G, is? 
E] = R] -¢] (4-73) 
u} = [co] [u] Ej (4-74) 
c] =[e.] [vw] м) (4-75) 


A phase plane portrait of the error and the derivative of the error 


u 


often aids an investigator in the analysis of a nonlinear system, By 
the use of the backward difference formulas of numerical methods, the 


Di 





derivative of the error may be obtained using the past values of the 
discrete error, The derivative of the error may be calculated by use 


of the present and one past value en 


zA , x 
е! = $ (e, = ê + O'T) (4-76) 


where the inaccuracy is of the order of magnitude of the sampling period 
T, If the present value and two past values of the error are to be used, 


one has °° 


e1 = (Зе = de.) + 4e,_5) + 0(12) (4-77) 


Use of three values of error has reduced the inaccuracy of the approxi- 
mation to within the order of the sampling period squared, In order for > 
this calculation to be accurate, the sampling period must be short with 
respect to the time constants of the system, This requirement is not 
түбітінен since this condition is necessary for stability in a sampled- 
data — system, 

If it was useful to postulate the phase space and determine the 
derivatives of higher orders it is possible to use the following formu- 


las, the choice of formula depending on the accuracy р 


// 
Br = = (е, == 2e 1 + Pio) * O(T) (4-78) 
m 
- = = (е, = Зе 7 + eo = en-3) + O(T) (4-79) 
77 1 

ог = з (2e. = ecl * 4e 2 == en-3) + O(T?) (4-80) 


Formulas for higher derivatives using backward differences are readily 
available,’ 
Writing the above formulas in a matrix form results in the follow- 


ing matrices with an accuracy of the order of the sample period, 
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1 0 0 
кз 0 0.. 
E - Gp 1 1 0 О E] (4-81) 
o o ere = 
| 1 O» . o 
Й КО 0 
E < (|1 -2 0 Е] (4-82) 
0 1 -2 1». 


The matrices with an accuracy of the order of the square of the sample 


period are; 


3 0 0 о о 9 
2 о о 
7 D 1 4 al 3 0 0 | 
|= Glo 4 4 0 o E] əн 
puc ms 5 o. 
2 0 0 ORE 
ES 0 0 
77 
Bl= (qn) 4 FF 0° 9 | дең) 
O =l 4 =5 za 


4-6, Analysis of Continuous Control Systems 
The analysis of continuous control systems can be accomplished with 
the use of time-domain infinite matrices, This is possible through the 


introduction of a fictitious (mathematical) sampler or samplers in the 
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continuous closed loop, The approach is basically that of the use of 
numerical analysis methods in the solution of a differential equation, 
Therefore, the accuracy of this approximation depends upon the sampling 
rate of the fictitious sampler, However, the fact that it is not neces- 
sary to calculate the closed-loop roots or to invert the closed-loop 
response equation in order to obtain the time response, is of great im- 
portance, 

Consider a continuous system to be investigated, which is a single 


loop system as shown in figure 4-20, 


м 


C (s) 





Figure 4-20, Continuous Feedback Control System 


The approximation introduced by the fictitious sampler will depend upon; 
1) The location of the fictitious sampler 
2) The form of the fictitious hold circuit 
3) The frequency of the sampling 
Since it is important for the sampler frequency to be many times greater 
than the highest frequency of the input signal to the sampler, the loca- 
tion of the sampler is usually chosen to be in the feedback loop before 
the H(;) as shown in figure 4-21, This location takes advantage of the 


filtering of the input signal R(s) by the plant G(s), 
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C(s) 





Hold 
Network 


Figure 4-21, Location of the Fictitious Sampler 


The fictitious hold network is present in order to reconstruct the 
continuous signal from the sampled signal. In actual sampled-data sys- 
tems a zero-order hold is usually used for its practical realizability 
and for stability considerations, However, for a fictitious hold the 
possibility of straight line and parabolic approximations should be con- 
sidered, Consider a straight line approximation hold which can be 


achieved by a triangular hold circuit which has a transfer function!” 


-sT,2 sT 
б, (s) i (1-е е 
old s?T 
The approximation of a time function is shown in figure 4-22, This. hold 


is not physically realizable, but is very useful for mathematical appro- 
ximations, 


ТШ) 





Figure 4-22, The Approximation of A Time Function 
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The sampling frequency determines the accuracy of the mathematical 
approximation and must be balanced by the amount of calculation that will 
be acceptable to the investigator, As the sampling period approaches 
zero, the inaccuracy is approaching zero, but the number of calculations 
rapidly approaches infinity, Fortunately, it has been determined experi- 
mentally that there exists a reasonable sampling rate for closed loop 
systems which yields less than 5%error in approximation, For a type I, 
second order system, this approximation holds when the sampling frequency 
is ten times the magnitude of the pole of the plant G(s), As an example, 
for a plant with a transfer function; 


K 


G(s) = ЫТ 

the sampling frequency should be 10 radians per second or the sampling 
period is approximately one-half second, The calculations of previous 
sections were carried out with a period of one second and a type one sys- 
tem, Therefore, the number of calculations necessary to determine the 
approximate response to a continuous-data system are not impractical, 
The simplest method of determining the necessary sampling rate is trial 
and error, One calculates a few points on the output response using a 
trial sampling rate, then recalculates these points using twice the 
sampling rate, If there is a neglible change in the results, then the 
former rate was sufficient for the desired accuracy, 

This method of approximation may be used on non-linear or time- 
varying systems and therefore all the previously discussed methods of 


analysis will apply. 








CHAPTER 5 


INTRODUCTION TO THE DESIGN OF CLOSED LOOP CONTROL SYSTEMS 


5-1, Introduction 

The method of time-domain infinite matrices may be used success- 
fully for design of closed loop control systems, The design or develop- 
ment of control systems for a specific application is a practical pro- 
blem of great importance, The designer may readily apply the analysis 
methods of the previous chapter to the design problem, They are used 
to evaluate the performance of the system under the specified operating 
conditions, Also, with the aid of experience in calculating the response 
directly in the time-domain, the designer may determine what system 
parameters must be adjusted and in what manner, 

Given the basic specifications of the system and the performance 
requirements, the designer often must determine a compromise between 
conflicting requirements, In the design of closed loop continuous sys- 
tems, many designers rely on the open and closed loop frequency response 
curves of the system to indicate the performance of the system, These 
frequency response techniques are not very useful in sampled-data sys- 
tems, Correlation of the time response of the system with the frequency 
response is not readily achieved since a transformation of the z varia- 
ble, into a new complex frequency variable w, is necessary to map the 
unit circle of the z-plane into the entire left half of the w plane, 
Therefore, the use of frequency response characteristics, such as the 
height and frequency of the resonant peak, and the bandwidth, is very 
limited, Design by use of the root locus method in the z-plane is 


limited since 1) the design only considers the response at the sampling 
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instants, and the intersampling response may be wholly unacceptable, 
2) Correlation theorems between the time domain and the z-plane are 
accurate only under certain conditions, 

The time domain performance indices such as rise time, settling 
time, over-shoot, and number of oscillations are readily applied to 
sampled-data systems using the time domain matrix method, Furthermore, 
this method allows one to design directly in the time domain and deter- 
mine the compensator necessary for the specific application, This method 
of design of a digital compensator shall be discussed in chapter 7, 

For the design of continuous systems, the designer has all the 
standard techniques at his disposal for the selection of a system adjust- 
ment or compensator, The methods of Bode, Nichol, and the use of the 
root locus all may be applied in the design of the continuous system, 
Then, the sampled data approximation may be introduced and the design 
evaluated directly in the time domain by means of the time domain matrix 
method, Therefore, all the background of previous determined theory and 
methods will apply profitably, The designer may use these techniques 
to determine the type and location of the compensator, and the parameter 
values, 

The designer may determine the necessary compensator directly in 
the time domain by use of the digital compensator design technique dis- 
cussed in Chapter 7, If the system is a continuous system approximated 
by a sampled-data system, the digital compensator may be converted to a 
continuous compensator by time domain network synthesis techniques, 

In addition to these design problems, one may consider the design 


of adaptive systems, Due to the great interest of these systems, they 
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are reserved for discussion in Chapter 6, This chapter shall be con- 
cerned with the use of standard design techniques and the use of time 
domain matrices to design various types of control systems such as non- 


linear and time varying systems, 


5-2, The Design of A System Utilizing Time Varying Gain Compensation, 
Consider the design of a closed-loop, unity feedback sampled-data 


system as shown in figure 5-1, The basic design steps are, 


Ше re) | Hola Plant e 
© G(s) 


Figure 5-1, Unity Feedback Sampled-Data Control System 


(1) Determine the response for the given system G(s) and compare 
with the required performance specifications, 

(2) If the performance is not satisfactory, adjust the gain or 
choose another plant if possible, Otherwise, introduce a compensating 
element or block in the closed-loop, Choose the compensating component 
on the basis of factors such as design criteria and design experience 
and the time response obtained for the uncompensated system, 

(3) Evaluate the compensated system response and readjust the sys- 
tem parameters if necessary, 

In this section it shall be assumed that 1) and 2) have been accom- 


plished and it has been decided to attempt to compensate using a time- 
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varying amplifier in the forward transmission path as shown in figure 
5-2, The variation of the gain must be chosen to provide the desired 
output response performance, Changing the gain of the amplifier does 
not alter the system dynamics and therefore the achievable results are 
limited, The designer may alter the response, but it is only possible 
to compromise between desired performance indexes, For example, it is 
possible to reduce the rise time, but only with a resulting increasing 


maximum overshoot, 


R (s) e E (s) 


Plant C (s) 


G(s) 





Figure 5-2, The Compensated Control System 


An example will best illustrate the possibilities, Consider a sys- 


tem with a transfer function of 
1 


= SD (5-1) 
and a zero order hold and sampling period of one second, The input sig- 


nal for this example will be a unit step, Then the uncompensated output 


at the sampling instants is found to be, 


G E С 
0 0 
3679 079 
„7575 1,00 
б] = | ‚45 1,40 
| ,9685 1.40 








Therefore, the uncompensated system has a rise time of two seconds and 
an overshoot of greater than 407 occurring between the third and fourth 
sample, If the specifications call for a maximum overshoot of 30% and a 
rise time less than four seconds, the required response can be achieved 
with a time-varying gain, The designer learns, from the experience of 
calculating the response directly in the time domain, that the first two 
error samples largely determine the magnitude of the maximun overshoot, 
Therefore, the designer choses an emplifier with a gain of one-half at 
the first two samples, and a gain of one thereafter, This amplifier can 
be practically realized by constructing an amplifier which switches to a 
gain of ,5 when a step input is applied, and switches to a gain of 1,0 
after two seconds ( n = 2). This system will give a more desirable res- 
ponse than simply lowering the gain to one-half for all time, This fact 
shall be verified later in this section, Calculating the compensated 
response one obtains, 

c]= [a] [a] E] = fc] F] (5-3) 


and therefore 


0 E О о о. 1 
„3679 0 ,5 0 0 „816% 
„7675 О 0 1 o. И 

С] = | .9145 of o o0 F, (5-4) 
9685 > y í 

| ¿9884 | | 
0 E 0 
„3679 „4080 1639 
„7675 4662 „5338 
= | „9145 „0581 | = „9419 (5-5) 

,9685 -.2265 1.2365 
9884| |-.2733 1.2733 
КОШЕ} | 1.1238 
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The intersample ~zsponse is found as follows, 


Eta .5 
| | .01 4080 wll 
c(m)| =[c(m)] P| = | ‚8590 4662| = .7307 (5-6) 
| 9481 0581 120607 
9809 -.2265 1.2884 
В .2930 | -.2733 1.2137 


It can be seen that the maximun overshoot is about 29% for the compensated 
system and the rise time is 3.15 seconds. The response of the uncompen- 
sated, compensated, and the system with a gain of one-half for all time 

are shown on figure 5-3, A comparison of the rise time and maximum over- 


shoot is given in table 5-1, 


System Rise Time Maximum Overshoot 
(seconds) (Регсеп%) 


"часы 


Uncompensated (Gain = 1) 
Time-varying Gain 
Gain = „5 for all time 





Further discussion of the application of the time domain matrix method 


to the design of control systems is presented in the next section, 


5-3. Application of the Time Domain Matrix Method to the Design of 
Various Types of Control Systens, 
In this section it is intended to show the great scope of possibi- 
lities of design of various types of control systems using the time-domain 


matrix method, 


A) Single Loop Linear Sampled-Data Control System 
When the preliminary analysis of a linear sampled data system reveals 
that the overall transient performance is inadequate, compensation techni= 


ques must be emplcyed in order to improve the system performance, For 


a single loop system, the simplest and most direct step is to change the 
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- Time Varying Gain 
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system gain, Usually, however, this adjustment alone is not sufficient 
to satisfy the design requirements, Therefore, it is necessary to in- 
sert a compensating network in the system in order to achieve the desired 
response, This compensating network may be cascaded in the forward or 
feedback channel or inserted as a minor feedback or feedforward loop. 

The sampling process in the control loop complicates the choice of the 
location for the compensating network, For error sampled systems, 1% 

is found advantageous to operate on the sampled error by a cascade com- 


pensator as shown in figure 5-4, 










R(s) E(s) 







Controlled 
System 
G(s) 


Compensating 
Network 
G, (s) 






H(s) A 


Figure 5-4. Continuous-Data Compensation 


This compensator operates on the sampled and held error 3nd yields con- 
tinuous data information to the controlled system, Furthermore, for the 


output response one may write, 


c] = [60,6] Е] (5-7) 
where | GGG] is the system matrix for G(s )c, (s)G (s). Investigation 
of this form of compensation reveals that it is difficult to stabilize 
a sampled-data control system containing higher-order integration with 
the use of linear continuous-data networks, The stabilization and coit= 


pensation of sampled-data systems by means of continuous cascade 
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compensation is further complicated by the calculations necessary in 
order to find the matrix for the overall system for every trial compen= 
sator, This is not a disadvantage peculiar to this method, but also 
results when using the z-transform, The use of a continuous-data net- 
work for compensation is not to be excluded, These networks are simple 
R-C networks coupled with amplifiers and are easy to realize, Further- 
more, for simple systems, they are sufficient and therefore probably 
desirable, However, for more complex systems it is often necessary to 
evaluate many trials in order to arrive at a reasonable compensation 
design, 

Therefore, it often becomes desirable to use a sampled-data network 
as a cascade compensator as shown in figure 5-5, Then, for the output 


response one obtains, 
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Figure 5-5. The Digital Compensator in A Sampled-Data System 


c] 


where P| 


[one] [o] £] = [ac] >) (5-8) 
L5] £j (5-9) 


Therefore, the sampled input to the controlled system is the output of 


the compensating device which is a result of a transformation of the 


error sequence, As was stated in Chapter 3, there is always a [D] which 
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will transform a diverging Sequence E] into a stable, converging se= 
quence P| „ Furthermore, the design of [D] is relatively easier thar 
the design of a continuous data compensator since the (6,6) matrix need 
be evaluated only once, The sampled-data compensating network can be 
realized as a program of a digital computer or more simply as a data- 
processing network, It has been shown, that a digital processing unit 
may be realized by operational amplifiers and electronic ns 5 
Therefore, a digital compensator may be part of a large computer program 
for missile control, or simply a controller for a DC motor, The realiza- 
bility and the synthesis of digital compensators shall be discusseú 


further in Chapter 7, 


В)  Multiloop Sampled-Data Control Systems, 

The design of multiloop sampled-data control systems is more com- 
plex and difficult than design of single loop control system, The use 
of linear continuous data compensation networks has the same limitations 
as discussed im the previous paragraph, Therefore, the sampled-data 
compensating network is used more often, A fundamental problem in the 
design of multiloop systems is the selection of the location of the com- 
pensators, This problem is difficult to solve in continuous systems as 
it is in sampled-data systems, One advantage of the time domain matrix 
method is the availability, in the calculations, of the response at 
intermediate points throughout the multiloops, wherever a sampler exists, 
Knowledge of each response allows the designer to use this information 
to adjust the parameters of the compensator, If there is only one 
sampler present in the multiloop system, then this advantage is not pre= 


sent, If a rulsed-data network is used as the compensator, then often 
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feedforward control wil) aid in the elimination of disturbance innuts, 
Consider figure 5-6 which shows a feedback control system with a disiur- 
bance input and feedforward control, Then the output due to the innut 
signal C, and that due to the disturbance signal Cu can be considereä 


separately as, 


چن سن 
U(s)‏ 
x‏ 






: l ls! 
= fm Hol“ 
G (s )=6, 6, (5) 
Figure 5-6. Control System with A Disturbance Input 
J =] (5-1с) 
Cul = (Ы + [D,] [c] uc] (5-11) 
where Р| = РИ R | + [ D, | R] 


Then, the output due to the reference may be written as, 


c| =[e] (To, Ф [o4] fi 1] + [D,] [ep R] (5-12) 


Therefore, the output due to the disturbance may be minimized by means 


of D,, and the digital network D, used to design the system for the out- 
put performance with respect to the signal input, Multiloop systems wii! 
be considered further in the next chapter, especially the conditional 


feecback system, 


C) Nonlinear Control Systems, 


The techniques used in analyzing nonlinear systems as presented in 
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the previous chapter are very useful in the design of nonlinear systens, 
The nonlinear component is considered as a time-varying amplifier in the 
system where the gain is dependent upon the input magnitude, The design 
of linear as well as nonlinear compensators for the nonlinear systems 

is possible, Also it is usually an advantage to have available the sig- 
nal magnitude at various intermediate points in the system, 

The first method of compensation to be considered would be the in- 
sertion of a linear compensating network at a point in the system located 
before the nonlinear element, If this proved to be unsatisfactory, then 
a pulsed data network could be inserted before the nonlinearity in the 
loop, Also, since the nonlinearity is considered as a magnitude sensi- 
tive device and treated as a time varying amplifier, an interesting possi- 
bility for compensation would be the use of a time-varying amplifier with 
essentially complementary characteristics to that of the nonlinear device, 
If the amplifier and the nonlinear element were cascaded, then the magni- 
tude of the signal into the nonlinear element could be kept within the 
linear region of the component, In the case where the nonlinear ET 
has storage of energy, this could not be achieved by a time-varying am- 
plifier, Therefore, in most cases a pulsed-data network is used as a 
compensator and designed on the basis of performance criteria and evalua- 
tion, 

The design of a nonlinear system usually involves some trial and 
error steps in order to evaluate effective design changes. The design 
of a relay servo, for example, may require a number of trials in order 
to arrive at the proper output voltage on the relay and acceptable values 
of dead-zone and hysteresis, The use of a nonlinear element as a com- 


pensator also may require many design trials in order to arrive at an 
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acceptable design for all the desired criteria, If there is a reason for 
using a nonlinear device as a comnensator, the calculations are as easily 
accomplished as for the inherently nonlinear system, Furthermore, the 
shape of the nonlinearity necessary to give the required output response 
nav be determined, There are no general restrictions imposed on the de= 
sign of nonlinear compensators, This is a distinct advantage over pre- 
sent methods of design of nonlinear compensators, The method may be 
applied to the solution of systems with multiple nonlinearities with the 


same eäse, 


D) Continuous Data Control Systen, 

The design of continuous data control systems may be achieved by 
the use of the time domain matrix by the introduction of the approxima= 
tion of the fictitious sampler and hold, The approximation and the atten- 
dant error is discussed in section 4-6, Basically, the sampling rate 
must be sufficient and a fictitious hold introduced in order to produce 
a continuous input to the controlled system. A continuous-data compen= 
sation network may be selected by any of the standard design techniques 
and then introduced into the fesáforward channel and its compensating 
effects evaluated by the time domain matrix method, A single lcop con- 
tinuous system with fictitious sampler and hold is shown in Figure 5-7, 
The fictitious sampler and hold are usually inserted in the feedback 
loop in order to take advantage of the filtering action of the forward 
channel on the input signal. The use of lead and lag networks as con- 
pensating filters may be investigated directly in the time domain by this 
method, 


If it proves desirable to investigate the use of a compensator in 
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Figure 5-7, Continuous Control System With Fictitious Sampling. 


the feedback channel, the introduction of a second fictitious sampler 
and hold will allow for compensation by a D matrix in the matrix equa= 
tions, Then the discrete compensator D found to be desirable may be 
approximated by a continuous network in the time domain, By techniques 
of approximation in the time domain a network can be synthesized to 


259 ; 
*” The time necessary 


yield a prescribed output for a prescribed input, 
for the calculation of the digital compensator and evaluation of a suita= 
ble continuous network may be considerably less than that for a design 


carried out for the continuous data system by standard s-plane or fre- 


quency response techniques, 


E) Time-Varying Sampled Data Control Systems, 

The design of time varying sampled-data control systems is possible 
with the use of the method of time domain matrices, The standard me= 
thods of analysis and design for non-time-varying systems, that is the 
z-transform, and frequency response methods, are not applicable, The 
time varying system is represented by the G(n,k) matrix and allows the 
general design procedures of the preceding paragraphs to be applied, 


Furthermore, the design of a continuous Лаа time-varying system is 
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possible through the introduction of a fictitious sampler and hold, 

This approach offers great possibilities to the designer, in the analy= 
sis and design of systems of which the dynamics vary with some variable, 
An important example is the high performance jet aircraft, where the 
dynamics of the aircraft vary with altitude and therefore time, A com- 
pensator may be introduced and a design evaluated, This compensator may 
be a continuous data or sampled data network, The limitations, advan- 
tages, and purposes of each tyne of compensator are essentially the same 
as those for the non-time varying systems discussed in the pieceding 
paragrarhs, 

If the variation of the time varying element is relatively large, 
it often is impossible to compensate with a non-time varying compensator, 
In this case it is necessary to use a time varying compensator such as 
a time varying amplifier, or perhaps a time varying network, The systems 
of this class are often termed adaptive systems, It is useful to consi- 
der systems with time varying dynamics (proles and zeros of the system 
transfer function), time varying gain, and time varying sampling rate, 
Research on these concepts has been carried үт The area of adaptive 
sampled-data systems shall be considered further in the next chapter, 

Finally, it is possible to use a time varying compensator in a non- 
time varying system, For example, it is possible to compensate a sys- 
tem for specified criteria by the use of a time varying amplifier, 


This possibility was illustrated by an example in the previous section, 


e 





CHAPTFR 6 


THE ANALYSIS AND DESIGN OF ADAPTIVE CONTROL SYSTEM 


6-1, Introduction 

An important class of control systems are those for which para- 
meter values of the system change as functions of some independent varis- 
ble during the period of operation, An important example is the change 
in the system dymamics of a supersonic aircraft with altitude, Another 
example is that of a chemical process where a parameter may change as 
a function of the ambient temperature, Usually, a fixed invariable com- 
pensator will be adequate over a restricted range of operating conditions, 
If the system is expected to operate outside this region, then some other 
form of compensation is necessary, The compensator often used, varies 
or adapts to the changed operating conditions, 

At the present time there is no widely accented fundamental defini- 
tion of an estive system, although several &re advanced in the litera- 
ume “BS À system which changes or adapts a parameter of the con- 
trolled system to drive the actual performance towards the desired per- 
formance shall be considered adaptive, The two basic segments of an 
adaptive system are 1) the identification of the dynamics of the con- 
trolled system, directly or by means of a related variable; 2) the 
generation of an appropriate actuating signal for the controlled system, 
Identification of the dynamics may be accomplished, for example, by 
measuring the impulse response, by measuring the response to white noise, 
or measuring a related variable such as the output to a known input 
Signal, The adaptive actuator may be a nonlinear, time-varying, or a 


digital device, For complex systems, it is common to use a digital 
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computer in which сазе the system must be treated as a sampled-data 


system, 


6-2, The Use of A Model in an Adaptive System, 

One philosophy of design of adaptive systems incorporates the use 
of a system model in the input channel as shown in figure 6-1, Тһе 
output of the model is then compared with the actual output and the 
error used to drive the system towards the desired output, Another form 
of a closed loop system using a model is shown in figure 6-2, This con- 
figuration has been called a conditional feedback system, the feedback 
of a control signal being conditional on the result of the error between 
the desired and actual system output, The necessary conditions for the 
desired performance are determined by the model and the Adaptive compu- 
ter provides the actuating signal, The model may be a physical simula- 
tion or analog of the process or a mathematical abstraction manifested 
as a set of equations stored in a computer, The Adaptive computer is 
often a special purpose digital computer, but may be a nonlinear or time 
varying controller, 

The application of time domain matrix methods to the analysis and 
design of adaptive systems follows the methods discussed in the previous 
chapters, The investigation of a sampled-data conditional feedback 
system shall illustrate the, basic aprroach, Consider a basic system 
where G, = 1, the adaptive comouter is a direct connection, and the 
feedback H(s) = 2, as shown in figure 6-3, The sampling period will be 
one second and the plant is a type I, 2m order system for which the 
system matrix is set forth in Appendix A, The goal of the design shall 


be to obtain a system output response for shifting system dynamics which 
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Figure 6-1, The Use of A Model to Shape the Input Signal 
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Figure 6-2, The Use of A Model in A Conditional Feedback System 
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Figure 6-3, A Sampled-Data Conditional Feedback System 


approaches the output response for the desired system, The desired res- 

ponse is present in the loop ás the output of the model, X(s). This 

desired response is compared with the actual and the difference used 

to drive the output towards the desired response, For example, if 

а = 1 for the unvaried system, then the output for a step input would be; 
o] = Colwnn 0, .3679, 1.00, 1.40, 1.40, 1.15, „#94, „798, ° > o} 

where (с) is the transposed column vector, This response is shown on 

figure 6-4 as curve number 1. Furthermore, if the varying system mani- 

fests its variation in a shift in the pole a, over the range a = „50 to 

а = 2,0, the simple single loop system output would vary in the limit 

as shown in figure 6-4 as curve 2 and 3, For the conditional feedback 


system of figure 6-3, one obtains the following matrix equations; 


Е] 
с] 


Then, one may write equation 6-1 as? 


R] + xj - 2c] (621) 


[c] El (6-2) 


i 





1 - 200 


1.3679 = 2, 
4000 - 25. 
E] - 2.40 - 204 (6-3) 
2.40 - 20, 
2,15 - 20, 
1,294 = 20 


For the system with the model in the limit at a = 2,0, one obtains, 


б Е 
ШІ 3 0 
.2838 1.226 „2838 
„4708 . 3626 „8127 
б] = | .4961 ,0478 = 1.1761 (6-4) 
„4995 -.1838 1,2919 
. 500 -.631 142085 
„500 = , 3064, 1.052 


e | e | 2 
This output response for the system with the model is shown as curve 
number 4 on figure 6-4 and the response at the limit for the pole a = ,50 
is shown as curve number 5, It can be seen, that the system with the 
model gives a response more closely approximating the desired than would 
the uncompensated system for either the limiting pole magnitude of a = „50 
ог а = 2,0, The simple conditional feedback system illustrated here, 
yields a definite advantage for varying parameter systems, If it was 
necessary to duplicate the desired response more perfectly, it would be 
necessary to use an adaptive controller in the feedback loop as shown in 


figure 6-2, Ina sampled data system this would usually be a digital 
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network or D block, The design of digital networks shall be discussed 


in the next chapter, 


6-3 An Adaptive Gain System, 
One possible method of adaption for a system is to set a parameter 
of the closed loop system dependent on a variable related to a perfor- 
mance index, Опе proposal advanced in the literature considers the varia- 
tion of a compensator pole based on a measurec figure of TETTE The 
calculation of the figure of merit, such as ITAZ (Integrated Product of 
Time and Absolute Error), is usually accomplished by a special purpose 
computer and consideration of the signals as sampled-data usually follows, 
One form of an adaptive system would use the magnitude of the 
sampled error to control the gain of the system, There is no signal 
storage involved in a simple gain adjustment system and therefore the 
compensation improvement is limited, However, this system illustrates 
the possibilities of this approach. Consider the system shown in figure 
6-5. The input shall be considered a step function, and the amplifier 
gain a function of the error. Then, for the equations of the system 


one obtains; 








R(s)= $ Controlled 
System 

G(s) = k 

5 (5+1 





Figure 6-5, Adaptive System With A Variable Gain 
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O O а, * 
€ 4 E | | (6-6) 
c] = Le] 7] (6-7) 


For a system where K - l and the amplifier gain is directly proportional 


to the magnitude of the error with a minimum gain of ,20, one obtains; 


ШӘ 0:00 o | 1 1 
0 0 O | .6321 „40 
00 3200 0 о | „0854 „0171 
Oo 0 |-.2277| -,0519 
Р|= 0 .2%20 0 0|-,3283|- -,1078 
О .120 0 |-,.312 -.0973 | (6-8) 
О .2430 | -,243 20957 
O „20| -,1684] -.0337 
-.1135| -.0227 
апа 
0 1 0 
298179 | 40 267 
INS 0171 „9147 
‚9145 | -,0519 1,22€ 
c| =|.9685 | -.ı078 | = 1.228 (6-9) 
9884 | =„0973 Т0 
„9957 | —„059 1.242 
„9984 | -,0337 1.168 
,9994 | -.0227 1a 
| | 1,076 


The response for the uncompensated system is curve 1 on figure 6-6, 
while the compensated response is curve 2, The compensated system has 
less overshoot, shorter settling time and also a somewhat greater rise 
time, For this type of Seon! the designer has the choice of the mini- 
mum gain and the function of error that controls the gain of the ampli- 


fier, If the pole of the plant is expected to change within the range 
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E 


а = 1 to a = „5, then the designer may evaluate the response at a = „5 
for the compensated system, The resronse for the compensated system 15 
shown on figure 6-6 as curve 3 anc the uncompensated response is shown 
as curve 4, It can be seen that, again, the compensated system will 
have a smaller overshoot. less settling time and approximately the same 
rise time as the uncompensated system, 

The change in the output response effected by this compensation 
scheme may not be satisfactory for many purposes, It is not the intent 
of this chapter to treat exhaustively adaptive systems, but rather to 
illustrate the application of the time domain matrix method to adaptive 
systems, Another powerful approach would involve relating the gain of 
the amplifier to the derivative of the error, The derivative of the 
error can be generated from the error sample magnitudes by the method 
of backward differences, The use of the error derivatives to control 
the gain or any other system parameter will yield a more desirable res= 
ponse than that using the error magnitude, A system which uses the 
error derivative to control the sampling rate is discussed in the next 


section, 


6-4. An Adaptive Samrling Frequency System 

sampled-data control systems usually have fixed sampling frequen= 
cies which must be set high enough to give satisfactory performance for 
all anticipated conditions, It is useful to reduce the sampling fre- 
quency whenever possible in order to extend component life and allow 
time sharing of the digital components, particularly the digital compu= 
ters, It is usually desired to have an efficient sampler, That is. 


over a given time interval, fewer samples are needed with the variable 


80 








е Ra чм cc л 1 
| ІЗ! КЕЛЕТ 
Өе, БЕ TERRE s Sri ta ae Er EHRE 
i | 


ЕЕЕ ЕЕ ЕЕЕ ЫЕНЕН НЕЕ Ше 
in SE ee ыы EEE 
EERE Ер ы р ы ЕН 

е ы р АНЕ 
ЕНІНЕН Нан ННІ 
ЕНЕ ИН 
ЕН Нн. 


| [Ep | i [xp LI ЛШ ө Шш | 


Be Е AAA A AE AA OE. a: 


ЕЕЕ ГЕ ы ЕДЕ A AHH 
4141 ——— =- fH gS ANA AMIA Д 
ыо иие mr | ep a На а гг UE вораи а а 
т ЕШ тиши Шаш ai ЕЕЕ БІРНШІ РГ ар а ы е jx 
ИЛЕ ВЕБЕ uU et ашытты 
шик | чараа eer TE Е | реја ЖЕ уш ee eee, ae | a Se 
ШИШИ ү т TT TITTI TTT TEV TTI TTT TET 
РБК це гаити 
E ааа PFMNINENMIENNEEN EN ыы 

IRE ү ال س کر ا ےل لا ا‎ ые ы Eee 
ШШ К Г атаа аа ауар о ЕЕ Pe Pe Es 
ea ааа an шан E 
Ea ALPES LL ICA a E ЖШ шш шш ШШЕ |] 
---------- т --а---- ТИЕГІ ЕРТЕН БН БЕНЕН 
mess] 1 Ku EL Li rid. DEMNM 
Ennis aisi ы 1 Ile} Tape) Adel 1 ale. Fal aaa ld et Ia T i 
eee ИЕ Г үү ауиминининликиин EEE же аш 
ШЕ NERRR 21111222. CALAMA 
aaa ИИТ т ыыы шы. "1 
е ниши AA RERUM ee "НЕ 
TE I T I I LO е | 
tante ا للا للل‎ О اال‎ 

IH ل ااا‎ ЕССИЕЛИИНЕЕЛЕГ 20 
AR c Ca ANTI a IA A ta a 
a | | | || cad oes 111 1 араа ЕЕ ББ 
E E ERR LEENA ELT LEL TTL ala alta шт тш гү 
ИЕШЕ КИНИНИ ыт Ат APRA AAPP 
о N | | 
Se ЫЕ) | Ыыы RE EN DR АЕ RR NR GV RR ER DR DD DT TR ШЕ Гә. ыш ЬЕ STS 
APODO LL. CAREER DA AMA MAMA A RA E A o. A AE A 
m 
Hf NL ЕСТ еар A 
---- SRA BECELA PESO: AN expe Me A ا ا ا‎ 


mu Qu mE QE IE Бе Кусага АСА 
HH EE-F-FHEHEFEFHEEESEISSEECE REL REEL deseen doi alae n Dad 


— ж — — 4 — —— | . سے‎ т 


AAA UL LLLLILILLLITÓEDSI IN CEIT E NBN NE: 
кшш ye) || | жмент а 
ШЕСІ, 11121: اا ا‎ HRS COLECTA 
ARSS TT LA A ја || Е Е оша ааа та 
аа ааа аа ајы | | Е раја: а Гарч. ча мај otr El a NE 
‚КЖ ЖЕШ дА Л А el ET] Pr td a ТЕТ ККА ТЛЕ А ER CIA 
ad E ESI lolo ci 

ГЕН аты EE 


-= | == t 


анаа ааа ЗЕЕ раја ре а 
TT pS fl 
ae لا للل‎ LS ON | p ITem 
E! el E AAA A A a 
ieee O lA j| E Uh ET Га Г ee ee NN 
Za лт + у E E E тты 
nl LL BS SNA 
NEP ELE Ee ET т шаасан 
e A A нан ышана 
параша yt T| Pt ee SI 
ыы соп AAA E aaa 


Compensated and Uncompensated System Response. 


Figure 6-6, 


m M س‎ pt а 
ИШИНИН EI I |} | | | st fe НЕН et eee ees 
A RA LLLI AA ARA AAA NE 
CRA SET He Hs 
Ca TIA PEO See SR RSE eee 
жашын Г ГҮКҮ ЕТТТ т ЕЕЕ ШЕ жс 
emmys |1, | LL TIL TTLLLLI TA gemi 
шм! e ТТТ поша ааваа ааа 
JECTS TE ЕЕ ы 2 итшттттиишшишиши C D 
aaa | AAA AH |i | Ts T P'o ШЕ eee ИЕ 
a утка -— Pp pre pr] Li rre NN 
не 1 -H te РІ л ша MN UNT 
ШШШ T т | a ee A 
Aaa EI LET ED 
ЕЕЕ ЕЕ ТЕШЕТ ТЕТЕ LAAL DL. 
_ tit tt EP Ip Е ЕТТТ н aaa e 
e ы Жу м т т SES 
Шиш р тр L ү E150 t О ЕН 
Ин ут ЕГЕТТЕ аа 
| TOHPUSEENEENMENENNNNENNNHEEBERNEBN р ы к ан UD D а 





frequency system than with a fixed frequency system while maintaining 
essentially the same response characteristics, A study of an adaptive 
system which varies the sampling frequency by measuring a system para= 
meter has been accomplished, ^" It was shown experimentally, that a samp- 
ler whose sampling period is controlled by the absolute value of the 
first derivative of the error signal will be a more efficient sampler 
than a fixed frequency sampler, 

Analytical methods of investigation were not available and in order 
to compare the experimental results with calculated results, it was neces- 
sary to develop a method of investigating sampled systems with varying 
sampling rates, Fortunately, tne time domain matrix method with a time 
varying system matrix may be extended to investigation of variable fre= 
quency sampling, Consider the system shown in figure 6-7, The sampling 
frequency is controlled by a function of the error, and could be con- 
trolled by a function of the sampled error if this was the available 


error signal as in a radar system, 
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Figure 6-7, Variable Samvling Rate System 


е2 





As was shown in the previous section, the derivatives of the error шау 
be obtained from the sampled error by backward difference formulas, The 


outnut response may be written as; 
сбх) | =| ai) | E) (6-10) 
where | 6 (n,x)] is the time varying system matrix 


and C(x) | and Е(к)| are time varying column matrices, That is, with 


the period of sampling changing with time, one cannot correctly write 
C | = (6,6, 02,0 4,04; e о о i (6-11) 


where the output samples are equally spaced T seconds apart, When tne 
sampling rate changes, samrles occur at various times and in this case 
k equals the number of seconds elapsed fromthe time origin, 

As an example, consider à system with two values of sampling pericd, 
T = 1 second and T = 2 seconds, Obviously, it is more efficient to 
have a sampling period of two seconds when the error is changing slowly, 
therefore the first derivative of the error may be used to switch the 
sampling Pe -iod from one second to two seconds, Consider a step input 


and a system with a transfer function. 


G(s) = Sn (6-12) 
Furthermore, for illustration, the sampling rate shall be considered 
to switch at the first sample when the derivative of the error is zero, 
which is at the first overshoot peak, Considering line A of table 6-1 
which is the response for this system with T = 1 second for all time, 
one can expect the same values of response for the first three seconds, 
At the output peak overshoot the sampling period switches to T = 2 
seconds and one obtains for equation 6-10; | 
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0 0 О е(0) 

g(1,0) 0 0 е(1) 

g(2,0) g(2,1) 0 е(2) жа. 
0 (%)| = ЛО) gl3,1) elB,21 0 e(3) (6-13) 

g(5,0) g(5,1) g(5,2) g(5,3) O e(5) 

g (7,0) g (7,1) g (7,2) 


, 
g(7,3) g(7,5) e(7) 
| e (9) | 
Therefore, using the values for the En given in appendix A for T = 1 
second up to n = 3 and then for T = 2 seconds one obtains the following, 


where the star indicates the T = 2 condition? 


0 0 10 
2679 0 ‚6321 
ШЕЙ5 30790 0 „000 

(4) = „9145 .7675 679 0 , 400% | (6-14) 
47884. ,9685 „9145 1,1353 4 -.1465' 
9984  .9957 .9884 1,9830" 1,1353 + ,2917° 


„9998 „9995 „9984 1,9842° 1.8830" 1.1353" 


Therefore, the output response is, 


End 


SOY = 0 

C(1) = „3679 

C(2) = 1.900 | 
с(к)] = C(3) = 1,400 (6- 

Ba, = 1270) 

CZ) = „7083 

C(9) = „9931 


— 


It can be seen that there will be fewer samples necessary for this varia= 
ble frequency and it is therefore more efficient, The response of the 
variable sampling frequency system is given in table 6-1, and it is essen] 


tially the same as the fixed frequency system with respect to overshoct, 


A p + 


E 


B Output-VariableT же] 1,00| 1.40 | 1 ds Иб 


Table 6-1, Response With Fixed and Variable Sampling Period, 
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rise time, and settling time, This subject is treated further in re- 
ference 10, where it was shown experimentally that a reduction in the 


number of samples of twenty-five to fifty percent may be accomplished, 


6-5, The Use of the Time-Domain Matrix for Analysis and Design of Adap- 
tive Systems, 

Adaptive systems may be treated as time-varying control systems in 
almost all cases, Therefore, continuous or sampled-data adaptive sys- 
tems may be profitably investigated by the use of time domain matrix 
methods, In fact, there may be no other method than can be used except 
solution by numerical methods ultilizing the digital computer, 

One of the necessary steps in the adaptive process is the identi- 
fication of the system, that is, the impulse response or the transfer 
function, With the use of time domain matrices it is possible to de- 
termine the pulse response when the system is at rest as outlined in 
chapter two, If the system is not at rest at any time, then the iden- 
tification of the dynamics must be accomplished by some other зе 82 

A predictor or learning system may be designed by the use of back= 
ward and forward differences, The evaluation of a forward difference 
equation allows the system to predict its next few values and adjust 
for optimum conditions, For example, the Gregory-Newton forward inter= 
polation formula may be programed in a digital network for the pre- 
dicted output as; 

с = 30 =- Зе с (6-26; 
using the present value and two past values requiring two storage ele- 
ments, The accuracy of prediction increases with an increase in the 


number of storage elements and if three past values were to be used 


eS 





in the calculation, then the formula to be programmed would be, 


cae 4C. - 6c, 1 + Ac» 2 (6-17) 


On the basis of the predicted value of the output response, the adap- 
tive system may adjust a system parameter or add a signal to drive the 
actual output towards the desired output, The accuracy of this method 
improves with an increasing sampling frequency as would be expected, 
This adaptive system using equation 6-16 and 6-17 may be implemented 
by the use of digital logic networks or operational amplifiers and el- 
ectronic NER: Equations 6-16 and 6-17 may be written in matrix 
form to aid in the calculation of the predicted values, 

The use of a special or general purpose computer in an adaptive 
system can be investigated in general by the method of time domain 
matrices, Опе adaptive system uses the integrated Product of Time 
and Absolute Error as the figure of merit and obtains the magnitude 
of this figure of merit by calculations on the sampled signals, These 
calculations are accomplished by a programmed numerical method and can 
be written in matrix form, Furthermore, the method of steepest des- 
cent may also be calculated by the matrix approach, In general, the 
use of numerical methods in the calculation of the figure of merit and 
adjustment necessary in adaptive systems may readily be accomplished 


through the use of time domain matrix methods, 
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CHARTER 7 


THE DESIGN Of DISCRETE-DATA COMPENS ATORS 


7-1, Introduction 

In the design of control systems it is desired to satisfy a set 
of specifications and response requirements, In many cases, this re- 
quires the introduction of a compensator in the control loop, The de- 
sign process may be carried out by the use of the time domain matrix 
for sampled-data systems and for continuous data systems by the intro- 
duction of the fictitious sampler, 

The introduction of a discrete-data compensator in the control 
loop will allow the designer to achieve a required output response, 
The compensator operates on the discrete-data input and yields a dis- 
crete-data output, This compensator may be realized by a special or 
„и digital computer, a logic network, or a circuit com- 
posed of operational amplifiers and electronic samplers, 

The process to be controlled is usually a continuous process and 
therefore it is necessary to investigate the response of the system 
between sampling instants while carrying out the design, Design by 
z plane or root locus methods does not include this possibility and 
often results in unsatisfactory intersample system response, The use 
of time domain cesign techniques has two basic advantages; 

1) the design is carried out directly in the time domain 

2) the intersample response is accounted for in the design, 

A sampled data system with a discrete data compensator is shown 
in figure 7-1, The D (s is the sampled data transfer function of 


the compensator, 


e" 








> а c—— ^ o— x MU Lm (е кекке 
-— € — O9 « o CM-—- cl ee 8 cm 


өн ж om eI A emm gem cm —HÓ — 
^ ce Ei u ж ж» жн), Е 



















De. a ---ж- Әт zn 





- c —— AA O o 


ж” ез MA. PUN >> «иь 
— ч - 1 нн gu -- aa! A 
«инна as: 4 1 ыыы “e эшш ж «жн», ж " 
LAA | © 4 
- - наш ж. Фе * 









۰ 
TH 15 5 = ә 





Figure 7-1, A Sampled-Data Syster with A Discrete-Data Compensator, 


7-2, The Fhysical Realizability of Discrete-Data Compensators, 
Tne digital compensator operates on the input signal and yields 
a transformed discrete tine sequence output, The z-transform equation 
for this operation may be written for figure 7-1 as, 
P(z) = р(2)Е (2) (7-1) 
Then, it can be seen that this operation may te written as a matrix 
equation as follows, 


Bills [n] E] (7-2) 


where 


when the compensator is time invariant, If the digital compensator 


is time varying, then one may write, 


Р} = [56.39 | E] 


where 


ag 


u - i 


_ = ама” 
nn ғ кен (же > 
о wu o _ A 
-—— © a = тү ү ү i 

> a ape o - — — — 
et ІП 


ж о: -—- o a 


u 64 Oa 


шт 





«90978 0 
айо) ай,1) 0 
[n(n,e)] = |a(z,0) a(2,1) а(2,2)... 


% 
It is necessary to determine what restrictions are imposed upon the D 
matrix by the requirement of physical realizability, 

А system is said to be physically realizatle, if the output signal 
of the system does not depend upon future information of the input sig- 
nal, In the matrix equation 7-1, this is expressed as the requirement 
that all the elements of the D matrix above the main diagonal be zero, 


The truth of this statement can be seen from the following matrix equa- 


tion. 


)7-3( 15 90 ا 


Expanding the second row (n = 1), one obtains: 

p = deo + dge, % 4,6, (7-4) 
Therefore, P, would depend upon the future input signal e, which is not 
available at the second sampling instant (n= 1%, Therefore, it is 
necessary for da = 0, and all elements above the main diagonal to be 
equal to zero, 


Furthermore, for a stable compensator, it is necessary for, 


d 
limit n+l 
п ә oO а, < 1 (758) 








If this limit is equal to one then the D block is a steady-state oscil- 


lator, 
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The discrete data compensators may be realized yt! 
1) Digital Frogramming 

2) Delay Line Networks 

3) Discrete-data RC networks 


4) Analog Computer нее 


7-3. The Sensitivity Matrix of the Sampled Data System with A Discrete 
Compensator, 
A logical quantative measure of the control property of a feedback 
system is the sensitivity which is defined as the relative change in 


the system transfer function T, divided by the relative change in the 


plant Р In the z transform notation, this may be expressed as, 
ا‎ _ ат(»)/Т(®) 
S #56) = ас (7-6) 


where me) = pe 


For figure 7-1, one obtains, 


6) = ЫТ (7-7) 


ы 2 
5 = 11607506 T 


This system possesses two degrees of frecdom, which permits indepen- 


and 


dent realization of T and S by means of the G and D components., Equa- 


tions 7-7 and 7-8 may be written in matrix form as, 


[7] = [в] DM [x] * (8] (91 7 (7-9) 
[s] =411] + [9] [o]] 7 


Then the sensitivity specifications fix the sensitivity matrix and 


and 


therefore, 
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fell] = (sI? -[x (7-10) 
The required output specifications fix the T matrix and since, 

[1] =[0] [9] [s] (7-11) 
one obtains; 

le] =|r] Is 2 [o] = (2-12) 
Therefore, equations 7-12 and 7-10 may Le used to obtain the required 
system matrices G and D, If the plant G is fixed bv power and load 
considerations, then a new comrensator, D,, must be introduced in the 
control loop, 

The sensitivity of an uncompensated and the compensated system is 
shown in figure 7-2, The uncompensated system is a type one, second 
order system with a zero order hold, The compensation reduces the over= 
shoot for a step input, to half the uncompensated value while reducing 
the rise time by half a sampling period, The magnitude of the sensiti- 
vity response of figure 7-2 is greatest during the first three sampling 
periods and settles out most rapidly for the compensated system, 

The use of the sensitivity matríx or response, aids the designer 
in understanding the effects of comrensation by a discrete data network, 
Furthermore, the variation of the system response with a system changing 
with time may be studied for adaptive systems by considering the sensi- 
ылуу, 

7-4. Design of a Closed Loop Discrete Compensator by Means of An Open 
Loop Discrete Compensator, 
For the single closed loop samrled-data system shown in figure 


7-3, a design of a discrete compensator is usually required in order 
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Figure 7-2, The Sensitivity of the Compensated and Uncompensated 
Я System | 
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7-3, Error Sampled System with A Discrete Compensator 


for the system to meet the necessary specifications, The output res- 


ponse at the sampling instants may be written in a time domain matrix 


=] [>] (r] «(el DTP s 2-19) 


This complex relation involving D causes the investigator to consider 


the design of an open loop compensator as shown in figure is 


Rís > 
(s) E 2 


7-4. An Open Loop Discrete Comrensator 





Then, the matrix equation for the sampled output resronse may be 


written as; 


с] = [e] [Po] R] (7-14) 


If an open loom compensator may be determined, then a closed loop com- 
pensator may be found, To find the relation between the open and closed 


loop compensator, equations 7-13 and 7-14 are equated and solved for D 
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yielding: 
[o] = [fir] - [»] [d] (7-15) 

It can be seen that this method takes the plant G into consideration іп 
the determination of the closed loop compensator, 

An example will illustrate the design procedure, Consider a sys- 
tem with a transfer function G(s) = star a zero order hold, 
and a sampling period of one second, А design of a compensator shall 
be accomplished for the best compromise for a step and ramp input, Тһе 


open loop compensator is chosen to yield a system output with a final 


value of one for a step input, Then, one has, 


зір ее 


(1 + 4,2 т К 





li 


Dg (z) 


2) (7-16) 


1 -1 - 
к (Landa dae 

This compensator will have a minimum complexity since the use of further 
delays such as а 27, increases the complexity of the closed loop com- 
pensator, The final value in the time domain, for a step input, ylelds 


a value of unity for the system output as expected, Тһе final value 


theorem is written and evaluated as follows; 


limit _ lim F 7 
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The output response for the sampling instants may be written for 


input as, 


E рак 


о 0 . 1 0 0 0 
Б 1 | 
с] = Sen) A ee A, (4-1) 1 0 
ee Бә 84 О -d, (4-1) n 
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Therefore, the output at the sampling instants is, 


En + a 
u € 5 


с = |. 
n 


п > 


The response between the sampling instants may Ре evaluated from, 


c] = [c] [2] в] 


Then, for half way between the sampling instants (m 


for the output response; 
e, (1/2) = 


For a ramp input, the output response may be written as, 
Bm ы [c] [»] R] 


R| = colum { Ot 2. ” 


Then, one obtains; 
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(7-18; 
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The output response is then easily evaluated for the step and ramp in- 
puts and the results for this example are shown on figures 7-5 and 7-6, 
The designer would be able to choose a value of d on the basis of the 
specifications for a step and ramp input, In this case, as a compromise, 
the designer might choose d, = -,60, Then, the matrix for the closed 
loop compensator D may be evaluated from equation 7-15, Since the sys- 
tem G(s) is second order, the closed loop compensator required will be 
third order, 

If more control over the output response is desired by the designer, 


0 
An 


then one may use a higher order Dg, For example, a higher order compen= 
sator might be; 


l 


-1 -2 =]. 
Do (z) = ca? (1*dz +4,2 5) 1 - ғ) 


(7-24) 
With this form of compensator, a deadbeat response to a step input is 
possible, Іп general, this method of compensation in the open loop, 
gives the designer strong control over the system output response, The 
method may be applied to linear sampled-data and continious systems, 


and time varying sampled-data and contirusus systems, 


7-5, Design of The Closed Loop Comnensator By Time Domain Evaluation, 
À simple method of the design of the closed loop discrete compen- 
sator should not be overlooked; that is, the arbitrary choice of a D(z) 
and an evaluation of the output response by the use of time domain ma- 
trices, The evaluation of the output response between the sampling 
instants avoids the possibility of hidden oscillations, The choice of 
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the D(z) may be aided by all the standard z-transform techniques, and 
evaluated by the time domain matrix, Tor the example, of the previous 
section, one might choose on the basis of z-plane techniques a compen= 


sator as, 


Ж Kiz = ,26782 
D(z) Т Z +d 


(7=25) 


Then, evaluating the outrut response for K= 2,4 = „50 апа а = ,72, 
as shown in figure 7-7, the designer would choose the compensator on 
the basis of the specifications, It is interesting to note that for 


K= 2, d = .72, one obtains 


E pM 2(2 - ,3679 
uo p SEET х а 


7258 


= de (7-26) 
This type of compensation in the z plane, commonly called cancellation 
compensation is misleading since the z-plane only accounts for the sampl- 
ing instants, The total resronse for all time may be evaluated using 


the time domain matrix. If Kis set equal to 2,718 and one obtains; 


G(2) D(z) = pay (7-27) 

and the closed loop, z plane response is, 
-— | 
Е) = ee Riz) (7-28) 


The response for this system is shown on figure 7~7 as curve number 3, 
It is obvious that it is necessary to consider the intersample response 


in most design problems, 


7-6. The General Characteristics of the Closed Loop Discrete Compensstor, 
It is worthwhile to look at the general form of the discrete comper- 


sator D, The compensator may be written in the z-transform as;> 
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D(z) = (2 + py) (z + poder = 


=] =2 
1 Ме + b,z > ооо 


К {1 Mri mm m .! (7-29) 


The time domain matrix is written as; 


28 O 
а, О 
mi= Kid, a T oO o (7-30) 
da da 4) 1 0 T 
| 
o p | 
о о | 
For a first order compensator one has; 
D(z) = HH (7-31) 


z2 +b 


where the zero is in the right half of the z plane, almost always the 


case, Then, one obtains; 
D(z) = K {1 - (a + b) l « b(a * b)z ^ - b^ (a « b)a 7? + ad (7-32) 


For this compensator, the matrix values are; 


а, = =(a+b) , d, = Ъ(а+ъ) , dy = -b” (a+b) y o...» (7-33) 


Therefore, if the values of d, and d, are obtained as necessary for com- 
pensation, the necessary D(z) may be realized as equation 7-31, If it 


is necessary to specify four elements of the D matrix then the compensa- 





tor must be second order, That is, if d,, d, d,, а, аге specified, then 


D(2) must be of the form? 


y. K (z=-a 2%с 
D(z) E z+b z+e) (7-24) 


The sampled impulse response of the discrete network weights the 


past and present values of the input signal, yielding the output signal 
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at the sampling instants, There are three possible samvled inpulse 
response time series for a stable discrete compensator, These possi- 
bilities are illustrated by figure 7-2 which also lists the necessary 
first order compensator and its z-plane pole and zero locations, An 
unstable discrete compensator occurs when b)1 in case c, where then 
the alternating series is a diverging time series, 

t is worthwhile to investigate further, the alternating weighting 
series of the compensator of case с, This discrete compensator would 
normally be used to stabilize the system while maintaining a desired 
rise and settling time, From equation 7-32 one can observe that if 
(a+b) is greater in magnitude than one, the second term of the con- 
verging series is greater than the constant term of one, Therefore, it 
is actually possible to weigh past samvled informat on with a greater 
value than present information, This condition may be used, although 
usually to weigh past data with the greater value would have an unsta- 
bilizing influence, A picture of the sampled impulse response of the 
discrete compensator is useful in designing a compensator as shown in 


the next section, 


7-7, Time Domain Design of À Discrete Compensator, 

The design methods discussed in the previous sections depend either 
on standard techniques on the z-vlane or on the choice of a trial com- 
pensator and evaluation in the time comain, These methods are powerful 
and are the ones used predominantly in practice, However, it would be 
of value to the designer if one could design the compensator directly 
in the time domain, Two methods of compensation design in the time do- 


main by the use of time domain matrices have been developed by the author 
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Figure 7-8, The Stable Sampled Impulse 
Discrete Compensator 
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and are presented in this section, 


А, Synthesis at the Sampling Instants, 
For the error compensated single loop sampled-data system shown in 


figure 7-3, one may write the following time domain matrix equations; 


awe Te] »] (7-35) 
P] = [o] E] (7-36) 
where E] = R] = | (7-37) 


Also, the equation Tor the response between the sampling instants may 


be written as; 
c| = [ew] r] (7-38) 


Therefore, if the designer specifies the values of the C matrix in equa- 
tion 7-35 on the basis of the required response, the values of the P 
matrix may be determined, Then, with a known input R, and a specified 
C matrix, the error matrix may be evaluated from equation 7-37, Since, 
P and E are then known, the required D matrix may be evaluated from 
equation 7-35, It must be noted that the designer has no control over 
the intersample response when using this method, When the P matrix is 
found on the basis of the desired output C at the sampling instants, 
the output between the sampling instants C(m) is then calculated from 
equation 7-38, 

Once the values of the sampled output are specified, the P matrix 
may be determined from equation 7-35, or by writing; 

ed (7-39) 

It is usually convenient to use equation 7-35 and avoid the inversion 


of the G matrix. It is possible to specify an output response at the 
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sampling instants which the compensated system will be unable to yield, 
in which case the calculated compensator, D(z), will be unrealizable, 
For most design problems, it is sufficient to find tre first few ele- 
ments of the D matrix and determine a D(z) comnensator on this basis 
from equation 7-32 or 7-34, 

An example will illustrate the procedure and the possibilities of 
compensation, Consider the system of figure 7-9 which has been considered 
in the previous sections, The output response of the uncompensated sys- 


tem is shown as curve number 1 in figure 7-10, 







R (s )=1/s 


T = 1 second 


Figure 7-9, Compensated Sampled Data Control System, 


It is desired to lower the rise time (for 90% of the final value) 
While lowering the maximum overshoot to 10 percent, Then, on this 


basis, set; 


T = column (о, DO des, Less as i (7-40) 
and using equation 7-35 evaluate the necessary values in the P matrix, 


For the first four sampling instants, one may write; 
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со О О О О Ро 
e, 3670 0 0 o | p 
с, | |.7675 .3679 0 O li De imn 
Gh | 9145 CTS HEN 0 P3 


=' = 


Using this matrix, the р, тау be evaluated step by step as follows, 


C4 


a E 
эы - XS «wies 


92-8 2Ро 1,0 -(,7675) (1,359 $ 
ИЕ. = ГЕ = “3679 = -,1170 (7-42) 
C3 -83Po -82P1 
o A lle 


Then, the D matrix may be evaluated on the basis of equation 7-36 and 


the calculated error values as follows; 


J = bo] 


1.35 8.0 


О 150 
- „1170 = К а, do 0 0 50 (7-43 ) 
UL аз dec .00 


Therefore, the d values are, 
Kdo = 1,359 9 Кд, = - „7965 , ка, = 2542 
and the D matrix may be written as; 
r О О 
к[р| = 1.359 | -.5e61 1 0 (7-44) 
%.1870 -,5861 1 
If a first order compensator is used to realize this compensator, equa- 
tion 7-31 through 7-33 will apply. Using equation 7-33, the compensator 


is founc from; 


а, 
(a+b) = ,5861 , b= = = есі „3191 


Therefore, a = .2670 and the compensator D(z) may be written as? 
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D(z) Е. 1,359 (2 -.2670) (7-4,5) 


z + ,2121 


1 =. 


M 


1.359 (1 -.5861z 5 « .le7027? - .05972 2 « „019127“ -•. 
For this compensator, the total response nay be evaluated using equations 
7-35 and 7-38, Тһе response of the uncompensated and compensated sys- 
tems are shown on figure 7-10, as curves number one and two respectively. 
Ihe output response between sampling instants is not directly taken in- 
to account in this design approach, Therefore, if too stringent require- 
ments are imposed on the response at the sampling instants, the inter- 
sample response may become unacceptable, The designer will usually eval = 
uate the intersample response for a design choice by using equation 

7-38, For example, if the specifications called for a maximum overshoot 
of five percent for the previous example, one might attempt to set c, 
equal to unity rather than 1,10 as was previously done, This design 


change will not effect the values of pg and p, or dg and d, that have 


been letermined, Therefore, 


ror > 17259 ытта о 
Са “бзГо “6214 
Ju uua ASE (7-46) 


£4 


Then, the gain and pole and zero of the first order compensator may be 


determined as follows; 


Kl, 2 pj -.5Kl, = -.4159 -.5(-.7965) - -.01765 (7-47) 


where K= 1,359 and K, = -.7965 as previously, 
Therefore, from equation 7-32; 

(а+ъ) = |а,| = .5°62 
апа 
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а, — „01299 


Б = ED = ee = -,02216 (7-48) 


6083 


bp 
il 


Then, the compensator may be written as; 


at = 6083) 
2 -.0222 (7-49) 


1,359 (1 ~ 58612 2 (01299272 - Оәданта”? - =» 4) 


D(z) 


The output halfway between the sampling instants may be calculated from 


equation 7-38 as; 


6 (1/2) | = I2) [o] Е| (7-50) 
neus 
6166 -.7965 

0/2] = [аә art Flo, 


ДЕТ - „0005 


mv Me „145 
_ 1.7531 | .5 |_ 5 
6744 | 0 1.0509 
es | 6 9305 | 


This response would have an overshoot of five percent between the sampl- 
ing instants and the design would marginally meet the specifications, 
This response is shown as curve number 3 on figure 7-10, ІҒ the desig- 
ner chose to set c, equal to 1,950, then the complete output response 

for a compensator calculated for this new set of response values at the 
sampling instants is shown as curve number 4 in figure 7-10, The maxi» 
mum overshoot for this compensation is 6,5 percent, A comparison of 

the rise time, maximum overshoot, and settling time for the uncompensated 
and compensated systems is presented in table 7-1, The designer would 
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chose the compensation that rrovided the closest approximation tc the 
original specifications, It can be seen that it is a powerful method 


of design of discrete compensators, The compensator may be located in 


Curve Nunber 1 2 5 

Maximum Overshoot 45 0,0 6.5 ж 

90 4 Rise Time de 1,65 1.65 seconds 
2.5 % Settling Time 16.0 IE 52 seconds 





Table 7-1, Indices of Response For the Compensated and Uncom= 
pensated Systems, 


the feedback channel or in one loop in a multiloop system and the design 
method loses none of its usefulness, However, if the designer has a 
set of stringent specifications, it may be necessary to account for the 


intersample response in the design, 


B, Synthesis of the Discrete Compensator Accounting for the Intersampie 

Response, 

It ís possible to account for the response between the sampling 
instants by the use of the time domain matrix equation for the inter- 
sample response which may be written as, 

c&)] = [ec] ]»] s] (7-51) 
Then, the designer may chose the first few elements of the D matrix by 
considering the sampled output while simultaneously considering the re- 
sultant intersample response value, Therefore, for example, the desig- 
ner will chose the Kig value of the р matrix by examining the resultant 
б. ааа с 5 values at the same time, This vrocess is carried on step by 
step until all necessary valuss of the compensator matrix are determined, 
In this manner, the next step would involve the determination of Ка, 


7 1 N © | choice Р ы ^ dae 
by a choice of c, and C1 6 The ch»ice of C and CLA? here á is 
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usually equal to one-half to yield the response midway between samples, 
is facilitated by the determination of an expression relating the two 
variables and its graphical presentation, The three basic matrix equa= 


tions are? 


d - [lo] =] =[0] fo] = о-ы) 
C (m)| = | Gta) | [o] E| = [2] | Gm) | Е | (7-22) 


g| » 8| «| (7-54) 


where order of matrix multiplication may be interchanged for time ir- 
variant systems, Then, for the typical system with gy = 0, and a unit 
step input, one obtains; 


És = Kd o &s (7-5%) 


апі E = Kigg 5 (7-56) 


It is usually necessary to choose Kd, such that the value of the first 
sampled response is less than the final value in order to achieve small 
maximum overshoot, Іп the previous example, for instance,  Kdg was 


set equal to 1,359, Then, one may solve for c 5 in terms of c,, 48; 





= B 5 6 е: 
с 5 572 C4 (7-77) 
4 Вл 


Equations 7-55 through 7-57 allow the designer to choose the value of 
Kdo while taking into account the intersample response, When the valve 
of Kd, is chosen, the value of e, is determined and then since 


do = 1 by the definition of equation 7-29, one obtains; 


K 0 1 
c] ee jE (7-58) 
Kd, | Ea 
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c(1/2)| = Р О (7-59) 
Кі, е5 e4 
Therefore, 
сә 7 Khg * K(g, * ge) (7-60) 
and 
e + EE ЕЕ СЕ u 


where the E,» К, апа е; аге known, Solving for c in terms of 2, 
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one obtains the general form of? 





g m 
с "di 5 р: + 9 (7-62) 
84 


where Q is a constant, In general, the linear relation may be written 


о 


ao 


1 





EJA 
) cr "9 (7-65) 


nen \ 89 n n 


It is often useful to show this linear relationship іп the form of a 
graph, 

Consider the example of the previous section A, where the specified 
maximum overshoot was five percent, The gain constant K shall again 
be set equal to 1,359, so that c, = „50 as in the previous example, 

The design procedure in section A did not account for the intersample 
response and it was found that this resulted in the responses presented 
in figure 7-10, The design of section A to satisfy the overshoot re= 
quirement would have a response as in curve number 4 on figure 7-10, 
Now, in order to carry out the design while accounting for the inter- 


sample response, one obtains the equations for c, and C7 e in terms of 
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1.359 | o 
с] = | ка, „3679 | .50 (7-64) 
ка, ‚7615 | e, 
апа 
1,359 | 2065 | 2 
c (1/2) | = | ka, m we (7-65) 
Ка, ‚es e, 


Therefore, one has, 
со 7 267958, + 142931 (7-66) 


апа 
су „ = .1065Ка, + ‚9104 (7-67) 


or 


Then сі 5 is related to c, as. 


С 


1.5 2 .2096с, ж ,5360 (7-68 ) 


and this relation is shown on figure 7-11, Considering the equations 
"7-66 through 7-68 and the figure 7-11, the designer might choose 


ca = „980 and с. = ,8234 as a compromise, Then, one finds Ка, = =,851, 


5 


and dy = -.626, The next step is to write the equation for c, and 


Cc, z Noting that the maximum overshoot for this design will occur in 


2345 


this interval, One obtains from equations 7-64 and 7-65; 


Cy = .3679KA, + 9648 (7-69) 

со g = 2106554, % 1,0194 (7-70) 
апа 

Co g T 289603 + «740 (7-71) 


Equation 7-71 is plotted on figure 7-11, and this curve shows that the 


overshoot may be limited to five percent, If c, is chosen as 1,040, 
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then c, = 1,041 and Kd, = +,2045, A compensator is then determined 


from equation 7-33 as; 


D(z) = 2ش‎ 6 (7-72) 


The total output response is then determined and is plotted in figure 
7-12 as curve 5. The uncompensated system response is again curve 
number 1, and curve 3 is the response for the design of section A 
which did not account for the intersample response, The rise time and 
overshoot are equal for each design, while the 2,5 7 settling time has 
been reduced 50 percent and there is no undershoot present, 

As another example, consider the single loop system as shown in 
figure 7-9, where the sampling period is again one second, the input 
is a unit step, and the system has a third order transfer function, 


The transfer function and zero order hold are written as: 


G(s) = AQi-e ) (7-73) 


s^ (s«1) (s«2) 
and then the system transfer matrix is given in appendix A, table A-1 
section II, Тһе uncompensated response of this system is unstable, 
aná the designer might try reducing the gain of the system in order to 
achieve stability, If the gain is reduced to two,or one-half of the 
unstable gain,the output response is shown in figure 7-13. This system 
has a maximum overshoot of 75 percent and a settling time of greater 
than 10 seconds, If the designer desired to maintain the same rise 
time while reducing the overshoot to 20 percert, it would be necessary 
to introduce a compensator, Using the discrete compensator design 


procedure, for a system gain of 4 as originally specified, one writes; 
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ШЕ. The Relation of the Intersample Response and the ` 
| Sampled Response 
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К 0 1,0 | 
M | Ka, 3362 | e, (7-74) 
Ка. 1,1869 ! е; 
odi, 9 
and x 
К ‚ose2 | 1,0 | 
o(1/2)|= | Ky | .1845 | e, (9-74) 
Ki. 1.4789 ег | 
Choosing K = 1,0 and therefore using the original system of 4, one 
obtains; 
C4 = .3362 
апа (7-70) 
e. = 663 


Then, one uses equations 7-74 and 7-75 to obtain the equations fer c 


C4 55 and the relation between them, as; 


с о 7 „33628, + 1,410 Kart, 
сі = „05228, + ,8221 "¡PR 
су; = 17320, + „5789 (7-79) 


Equation 7-79 is plotted on figure 7-14, Now, choosing c, approximateıy 


equal to 1,0, one might use d, 7 -1,U, Then, obtaining the relations 


for c4 and Co g One has, 
ез = «23624. » 1,0266 (7-80) 
655 = „05824, + 1,172 (7-82) 
= „зда F22 (7-82) 


Equation 7-82 is plotted on figure 7-14 and from this figure one may 


choose cy approximately equal to с, 5 &t about 17 percent oversnoot, 


c. 


7 








ho о, ж ө | ать а ee 


Б 
aa 
= _ а 


Then it is found that d, may be set at „5, A first order compensator 


is determined for these values and one obtains; 


| 2+ ‚5 elo. .-2 = - Y] 
D(z) o ena = (1 = © 22 = 2208 3 D „125% ^ де o) (7-841 


The response for the compensated system is shown in figure 7-13, The 
maxímum overshoot is 20 percent, the rise time is reduced to 1,8 secernis 
and the settling time is reduced to 6,4 seconds, 

These design methods are equally applicable for other system input 
signals, This will be shown by designing a compensator for the systen 
considered in this section subjected to a unit ramp signal input, Again, 
the sampling period is one second and the system and hold transfer fur 
tion i$. 


=з 
Ge) = | ee) (7-84) 


5^ (s + 1) 


The input signal and the uncompensated system response is shown on 
figure 7-14, The first two samples have a zero value, and one may 


write for the output response; 


On dies 
1 ЕК |59!) UD 
= کے‎ 
c | Ka, | .7675 | e, ЕД 


Kd, | „9145 | & 4 


| 
1 


and = 
K „1.065 0 


0(1/2)| = | қа, „61% | 1,0 (7-86) 
ка, 828 е, 


Then, choosing K = 2 so that c, = „735€, one may write for tne next 


sample interval, 


He 











х0! SH 


EUER 














m LLL jJ 
4 ЕНІН НЕ 
дее HHAH لان ر لل‎ 7 
БЕН ЕНШІ 


ее ШЕ а аш n 
„е Г ее е РЕ” 
Ы HH 




























Sie e TL 
pt DL III Ж ГЫ Шш ШШ ЕЙ А к 





"ГІНЕН. Т ПЕТ ГІТ ТТТ ТЕ x 119271 
Ela dale ANTONI TOI IE CCC oo аја E RR e iR 
ELA ABE CEA AA po ар Заа ашса ааа ac LE 







Reduced Gain [K = 2] 
| 
E. 
3 
E 
Е 
E 
Nu 
E 
[d 
[3 
n 
"m 
7 
E 
=== 
E 
1% 
MIL: 
a 
М 

| mr Jp т 
mi E Ұ 
| 
ү? 
Io 
E 
w 
m 
E 
| 
ы 
| 
| 
[u 
E 
Е 
Е 
E 
2 
E 
| | 
m 
E 
A 
Ж 
Ж 
a 
Ж 
E 
X 
x 
E 
b 
Ed 
a 
¥ 














ЕЕ | 

———— 
е | | | dL A 
Ер 
Кт E 
— – 













| 
! 
| 
| 










EX Е 
E 
HL 


meee cece 
سا للل ل الالام‎ 
ЕЦЕ 
nn 6 de 
eet tt REN NC EUER REIR 

| B 
Tee TTT 
ee 







EESEEEEENNUENNRENNEERNNEENNENANNNNENNNNER EM 
I-II LL qe AAA AA а 
Y ККК ТТЛДЕ ТИЛЕТ ы! 
ساسا“ ااا‎ Pa ET ied, RS Т А ЛЛ ШШ 
-2---------------- -- 


һы ! 
CA AR ЛУК АТ ү ғала 
A a ЫЫ МЫ ше аад ШЕ АД 
Ea d ol dd i ic) [mis te ae ala Fl dl e AI 


ЕЕЕ ЕЕЕ ЕЕЕ 
ЪЪ мыры ақылы ыы ҚЫ ЫЗЫ АА ы de ale le ale dle el e SS fet 
- ШАЛ мад) ЕНЕ На ааа o ЕЕЕ. 
Es ГГ ГГ!!! Г ГГЦ yS 
a ere te te К К ыт е у тд | al LA мий Be een DELI m 
ШИЕ а-а-а”. 
A AAA m 
4 ET AES AL NABERA AT FEE CAE A ELECTRA dan АЛҒА МЕЛ АТА. ІЛ алалы лара 
a е ыыы нын 
| | | EE HE 






= 










ы ыыы e 








--4 -- 


Compensated 























> = =r 










































E 
L- AA ы АЕ Т LET SE ET Eh а Б дс TE ES 

КЕТТ A Ere rer rss 
eee eee eee PTTL PEA LT ee Se A A 


I BR ТТЕ 1 IEA МЕР WEN UTE 
ЕН ааа Ны Dei 


El 


H | 

Em | | | [ | шаша ыш Ес ЕГ Ce pala aa e 
ТЕЕ PS ЕН км 
ET TILL LEI IT Ro I TIECELLEELLCEF TFT FF BER TEE ts 82144142 | 
HH HH HR а | 
-L сат di e e de ble dd dal AAA PA 


-d 








































Em 121111. КЧ ГА 
E аты COSER 
БЕ o LLL E EEL bancada э 


НЕН 





Ggim rai ә Res SCOT Hea 










à 


apensated and Redu “ac. 


ғ 
% Ұз 
а? 


— A 


2g! 


E 
































































IIIIII] = EERE BS ЕЛЫ ӨЕ ВИ БАИН ИИ ЕШ ШЕШ И Н 
к m CSI T1 ШЕСІ” ae PEE ere E ЕЕ ВЕ 
ЕЕ Наннан 

ETE БЕНЕН Ama 









A ы E 










































ee |||, EE a mc „| 
2 > 
a aa arpa PAT EA A cis) Stee 
BIEN шыл іші, 8 аша: з= 
e [al La A | a e aa ET Sr a mm c NM M E 




















LEI TT INI | 
HH 





S la A 0 m BES 
ГЕНГЕ ЕНІП SS 
— — 002: 
سل‎ E A ee ыт рр тар г арааг леш ee 

ЕЕЕ EIE TAHA RHETT GS 
E E E a a LET EIN TAT el METIA 
ETT TTT STS С, ытаа p ab dass E E AAA 

n 






















mime | Ll ee ee (EPR aMSN Ж FE) | ---- 
Ы- LE O ыы ші oa Mell WELL m 
| ИШИ ШЕШ ет т ЫЕ ALMAS TP Е. 
ВИ EEE E EEE ETE A ale Т e 7. IE nu TES 1 
Е Еш 
A (o MM 4 ee ----- A Шаты 
TTS ELF oa Сү Т ЫЕ Т 
кы ET EEE REECE EEE аа аа EE 
met tT |] A le leida! Lale Tesla! BA 


Mea ada der Ph hcdddo RAR RR RRR RR RR HTT 


| Figure 7-14, The Relation of. Intersample Response and s 
| | The Sampled Response 


120 





Je 
E 
۱ 
۱ 
] 
| 
| 
| 1 
| » f 
I | 
Е; 1 
I Т 
I 
| "! 
‘i 1E 
:% 
E 


ШІ” 


„eu >] | 





ү! o 
ЬУ І 
IM tit 








с, = ,3679Kd, + 2,465 (7-87) 


с 5 = 1065, + 1,5026 (7-88) 


апа 


"2.6 = 2289; > ,7977 (7-89) 


From these equations aml the nature of the response, one might set 

Ki. — -,50, Then repeating these steps fer the next interval, one 
might let Ki, = + „20, The value of id, chesen wouid depend largely 
on the overshoot limitation, In this case, there would be no overshoot 


resent, The total response is shown in figure 7-14 as curve number 
p iig 


қ 
с“ 
г 


two, The steady state error of the compensated system is sixty percent 
of the steady state error of the uncomrensated system, The compensator 


to yield this response is found to be? 


D(z) = 2122-00) (7-90) 


if it is required to have a shorter rise time, while permitting an 
overshoot less than five percent, it is necessary to raise the value 
of K to three, Then if tte usual steps are carried out, one obtains; 
K=3 , Kay =-1.0 , M, v,5 
Then, the total response is found to be as curve number three in figure 
7-15, The rise time E greatly reduced, while the maximum overshoot is 
less than two percent, Also, the steady state bens reduced to 42 
percent of that for the uncompensated system, The compensator reces- 
sary to realize this response is found to be; 
D(z) E. 2 167 
(2 + „50) 
If the input signal was expected to alternate between a ramp and step 
signal, it would be necessary tọ design a compromise compensator with 


3 
L. 
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Ay 








respect to the design requirements, This is песеззагу, since a com 
pensator te be cptinum with respett “е с дітет set cf conditicns, can- 
not be expected to be optimum Tor 2 thanged set of sperating conditions 
and requirements, 

It can be seen that the design sprroaches presented in this sec- 
tion allows the designer a wide latitude of the choice of specifica- 
tions, ana affords direct control over toe output. response, Furthermore, 
the time domain approach gives the designer a direct picture of the 
time response, therefore avciding inaccurate and unvieldy correlation 
theorems for the z plane, The desigr. method which accounts for the 
intersample response allows the designer te wield control over the tote. 
time response, It is obvious thet this iesign procedure may be pro- 
grammed in a digital computer, the desiguer supplying the specifica- 
tions and system transfer matrix as tbe izput to the computer, 

A discrete compensator may be üetermines for a continuous system 
which is approximated by a fictitcus sampler and held as discusseä 
previously, Tren the calculated D(z) may be synthesized by & continu- 
ous compensator in the time donata,? Thereicre, this technique is 
not limited to sampied-Juta systems, bus may be applied to any type 
of continuous data systen, 

The location of tne Ciscrete ccompensator AS not limited tc the 
error channel and way he located in any arbitrary loop in the system 
as feedback compensation , 

Therefore, two methods cf design, directly ir the time domain, 
have been presented, These design procedures are entirely flexible 


in arpplivation, accurate in calculation, and rapid in sclution, The 














design in the time domain gives the designer a complete insight into 
the total time response of the systen, No other existing design pro- 
cedure can be carried out directly in the time domain, nor can any 

existing design procedure provide the flexibility and accuracy of the 


time domain matrix method , 





CHAPTER 5 
CON? LUS TONS 
?-]. Summary of Results 

The aim of this dissertation was to present a new method of 
engineering analysis and design for complex control systems, This 
method is the time domain infinite matrix method, The formulation 
of the infinite matrix follows from the convolution summation of 
sampled data systems, The mathematical basis of the time domain 
matrix formulation is presented in a discussion of the applicable 
concepts of infinite matrices and sequence spaces, This method of 
analysis and design is applicable to both continuous data and sampled 
data systems, For continuous systems it is necessary to introduce 
a fictitious sampler and hold of sufficient sampling rate to effect 
an accurate approximation, 

It is possible to analyze and design linear, nonlinear, and 
time varying systemof the continuous or sampled data class, Sam- 
pled data, time varying systems may not be investigated by any other 
existing method, Furthermore the investigation of nonlinear systems 
is greatly simplified by the time domain approach, Multiloop sys- 
tems may be treated with ease and the signals at intermediate points 
throughout the loops are readily available, Also, systems with 
multiple nonlinearities may be investigated, for which there is not 
a presently available method of analysis and design, 

Two methods of design of a discrete compensator for a sampled 
data system are presented, These methods are accomplished directly 


in the time domain and allow for a compromise of specifications in 
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the time domain, Also the response between sampling instants is 
accounted for in one of the two design procedures, 

The time domain matrix method mar be readily programmed on a 
digital computer and therefore provides а rapid analysis and design 


technique, 


8-2, Further Conclusions 

It is an important advantage that the design and analysis of 
systems using the time domain matrix method takes place directly in 
the time domain and not in any transformed complex variable domain, 
This advantage aids the designer in widerstanding and controlling 
the time response of the system under study, 

The availability of the intersample response is also an advan- 
tage to the designer, so that a design may be accomplished which 
accounts for the total time response ard not that solely of the sam- 
pling instants, Therefore, the time domain matrix method has strong 
advantages over the commonly used z-iransform on both of these points, 

Since adaptive systems may be treated as time varying systems, 
the time domain matrix method may be applied in general, Further- 
more, a learning control system with a predictor may be realized 
and investigated by the use of time domain matrices, The investiga- 
tion of conditional control systems incorporating a model was also 
accomplished, 

For nonlinear systems, the time domain matrices provide a useful 
analysis and design technique, The characteristic of a deliberately 
introduced nonlinear compensator may be determined and investigated; 


that is, it may be specified by the design approach, Furthermore, 
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the derivatives of the error or any other desired signal may be 
evaluated readily by matrix metnoós and used to determine the phase 
space response to aid the designer in the investigation, The use of 
a fictitous sampler and hold in a c:n*znucus data system permits 
these techniques of investigation and design to be applied, 

The design of the discrete compersator by means of the time 
domain matrix method provides the designer direct control of the 
time response of the control system, Therefore, the designer may 
specify the required system response and determine directly, the 
necessary discrete compensator, The i=2lizability, stability, and 
sensitivity of the discrete comvpensator was investigated, to provide 
the designer with an insight into the general characteristics of the 
discrete compensator, 

In comparison with other methods cf investigation and design, 
the time domain matrix method is more accurate due to its numerical 
formulation, Furthermore, the accuracy and ease of calculation for 
the matrix method remains the same regardless of the order of the 
system transfer function, while the &couracy decreases and the dif- 
ficulty of calculation increases for other methods, such as the z= 
transform, The time necessary for calculation of the time response 
of a control system is considerably less than that of standard com- 
plex frequency techniques, Actually. ths time response may be de- 
termined as rapidly as the determination of instability or stability 
by standard techniques such as Routn!s method, 

The transfer matrix of the system may be easily determined by 


experimental pulse techniques tc an accuracy of at least three 
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percent. This method of system charasterization may be used for the 
identification process in adaptive systens, 

Therefore, it has been shown that the time domain matrix method 
is applicable to a wide range of control system problems, Techniques 
have been developed which permít the solution of problems which are 
not solveable by any other method, Examples of such problems are the 
analysis and design of systems containing more than one nonlinear 
element, and the analysis and design of time varying systems, for con- 
tinuous or sampled data, 

A method has been developed which permits the design of a com- 
pensator for continuous or sampled data systems, without the use of 
trial and error methods, It also has been determined, that the 
difficulty of the application of this method and the labor involved, 
is considerably less than for other known methods when applied to 
systems of reasonable complexity, The advantages inherent in this 
method are sufficiently great, that it should find wide application 


in the engineering analysis and design of systems, 


8-3, Future Work 

The time domain matrix method may be applied to the comlete 
spectrum of control problems, Therefore, the possibilities for future 
investigation utilizing this method are boundless, 

For all classes of control systems, it would be valuable to in- 
vestigate the problem of the location of the discrete compensator in 
the control loops, and determine and classify the relocation of the 
compensator on the time performance indices, 


For nonlinear systems, the investigation of the use of the error 
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derivatives and the signals available at intermediate points through- 
out the loops, would be very worthwhile, 

For time varying and adaptive system, this method would bring a 
new insight into the very difficult problems, It is possible to use 
the time domain matrix approach to investigate the method of steepest 


descent for adaptive systems, as one example. 
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APPENDIX A 


THE SYSTEM TRANSFER MATRIX 


A-1, Determination of the System Transfer Matrix Using the Z-Transform, 
As discussed in chapter two, the actos of the system 

transfer matrix G is an important first step in the use of the time 

domain matrix method, It is necessary to evaluate the values of Eps 

the values of the impulse response at the sampling instants, lf the 

system transfer function is given in the Lsrlace variable, it is pos- 

sible, by standard methods, to transform to the z complex variable 

where 2 = E Then the values of the impulse response are obtained 

Бу inverting the z-equation to the time domain, The simplest method 


to accomplish this is use of divisicn to yield the response at the 


sampling instants, That is, G(z) may be expanded as: 


G(z) = Ео + Р4 2 + 552 P оосс ma 


(а-1) 
This series is the constant term and principal part of the 
Laurent expansion of G(z), . 


As an example, consider a system as shown in Figure A-1, 






Zero Order Hold 


let 
S 






R(s) ж 











G, (s) = 





T= 1 second 


Figure A-1, Open Loop Sampled Svstem 
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The z transform of this system may be written as, 

S], 
з к(1 ce ; ) (auo) 
© s?(s + 1) ) 


si 1 d m 
1% O E n 6792 + ,26422 “) 


m Ра 


G G (2) 


li 


Then dividing we have 


GG (2) 
K 


IN 


NS -3 
.367902 ^ + 676752 + „91453 7 + oes (а-3) 





This method is tedious in that cne must transform into the z variable 
and then divide, It also involves an inaccuracy introduced by the 
round-off in division, Jury” presents a method of inversion avoiding 
division which uses G(z) in the forn? 

=l =2 | _=N 
Po + р, 2 ғ рә? "^ 000 ра 5 


біз) 3 55-2 Ж. (a-4) 


+ + dle oo 2 
1 +4277 + 4272 +... 9, 


where m > me. 


Then, it can be shown that: 


бо = Po 
81 = р, 7 0480 (a-5) 
Бо = Pe = 04Е; = GOES 


апа 


Jury's method reduces the inaccuracies introduced by the division 
process, but remains inaccurate for a vaiue of n greater than four 
due to the roundoff error of multiplication, For a higher order sys- 
tem, the use of the z-transform aná inversion is unwieldly and inac- 


curate. 
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A-2, The use of the Impulse Response to Evaluate the System Transfer 
Matrix, 
A simple and more direct method of evaluation of the system 
transfer matrix is the use of the impulse response of the system, 
If the impulse response g(t) is determined, then the elements of 


the system matrix are found by substituting t = nT, That is, 


-— RAM, ; = g{nT) (а-6) 
=n 


Since, invariably the impulse response is made up of a sum of 
time functions, the arithmetic operation inyolved in the evaluation 
is addition instead of division cr multiplication, Therefore, the 
roundoff error of calculation is reduced as well as the difficulty 


of manipulation, Now recornsícer the example of the previou: section, 


o 


The system transfer function is? 


1 ¿ST 


-7 
Fea) ре 


GG (s) 


and the impulse response is; 


-l ( =1 -sT 
= j WA a e 
= La | L jatar} (a-8) 


S 
The p simply implies a delay cf cne sampling period, Then 


the inverse is. 


g(t) = a - (1 - | ^ г an) m $T 2 u(t-T) 


- 


where T = 1 second as on figure А-1, (a-9) 


Then the values of interest at the sampling instants when t = nT =n, 


are. 


g = T+ ro OUT - Tj 


-n ES - 1) 


1 +е E (a-10) 
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Then, it is found; 


Eo О 
Bo = NY = 226788 (а-11) 
go wm lt as - om = „76746, etc. 


Therefore, the values of the system matrix are rapidly and 
accurately evaluated using a mathematics table, This calculation is 


a great improvement in accuracy ever that of division in equation a-3, 


A-3, Evaluation of the System Transfer Matrix for Intersample Response 
Using the lodified-z-Transform, 

In order to determine the response of a system between the sam- 
pling instants, it is necessary to obtain the system transfer matrix 
for values of time between the samples, In order to accomplish this 
with standard notation, let m = 1 = A, where A is the percentage de- 


lay from the sampling instants as shown in figure A-2, 


g 
En+A an 


0 
> 


8 nel „AT و‎ 


а ао 5 


E rer 


Т ааа 
(n-1) i n+1 


Figure A-2, Delayed Sampling 
Then, writing the transformation t = (n =- A)T one has: 
t= (m = A)T= nel ro O<m< 1 


and the modified z transform is; 


p 
„э 
uh 








COD = ж. > gkT e anys (a-12) 


к= 


e 


The output response of a system between the samples Y(z,m) may then 
be determined for an imput X(z) as; 

Ү(2,п) = G(z,m) X(2) (a-13) 
G(z,m) may be determined and inverted by division or an alternate 
method to obtain the elements cf tne system matrix for intersampling 


response, Then the intersample output response is, 


Y (m) | = [cm | x | (a-14) 
where : 
£o (m) 0 Q ооо 
a(n] =| em) gom) o 
82 (m) g4 ќт) £g (m) ооо 
L o о о аж 4 
Since G(z) = 2G(z,m) (a-15) 
m—» O 


the | б] matrix may be determined by allowing m— 0. 

As an example of the determination of the modified system matrix 
using the modified z-transform consider the previous example, Тһе 
modified z transform of the system is figure A-2 is; 

G(z,m) = Канела s pog 4.(.2879ne = 7358) (2-16) 
27 - 1.90792 * ,36779 
where p, = 2.3679 = m(1,?679) a 


Expansion by the method of section A=1 yields 


Y" 


£o (m) = Ро (т) - m + e - 1 

gy (m) = p, m) = po (m)q, (а-17) 
bn г 2 

g 2 (m) = pa (n) = Py imag 8: rs (m) (a5 = qı ) 


and so on, 
1% 





A-4, Evaluation of the System Transfer Matrix for Intersample Res- 
ponse Using the Impulse Response, 


The values of the modified matrix are simple to evaluate by the 


impulse response method, One has, 


g(nT + AT) = gat) | (a-18) 


g, (m) 
lt = nT + AT 


g (nT + T = mT) 


Then, for the example previously discussed one has, 


(t =~] + "us 


Eo (m) 
t =ıT + AT 9 Те- 1 


Аа ыы 


(а-19) 


л ea А s Qu an 


Then, for example when A = 1/4 , m = 3/4 спе has? 
80 (3/4/- 1⁄4 = 3/4, = (7248 


For n? l one has 


ot -(t-T! 
~ a 


il 


е, (т) Т+е 


t = nT + AT 


+ е7 nts) > 67 0224-1) (а-20) 


Therefore, 


g, (m) = 1 + ur) rn 


к (ақа (а-21) 
кет iia Ч 


1 
H 


82 (m) = 


and form = 1/2 , л = 1/2 


с, (1/2) = ,61660 


1 
GQ 
ps 


(а-22) 


il 


ШЕ (1/2) 2595 


N 
29 





14. 





se. r—— — 

„ = Poe 
"218564 
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where A is another form of notaticn, 


А-5, The System Matrix for A Type I Second Order System with A Hold, 
The system transfer function for a type I second order system 


with a hold may be written as. 


=sT 
G G(s) = Ka =e (а-23) 


s^ (s * a) 


Then using the inverse Laplace transform the impulse response is, 





g(t) = > Е -(1- ty forts T (а-24) 
a 
and E 
e(t)= E (at = D 4 (aten E aem) tT 
a 
oe art nn] = 
a 
Therefore; 
E, = ЗЕ + 2 ЗЕ athe) | novi (а-26) 
а? 
апа | 
зе) = ا‎ a 2T (n+4) = | (4-27) 
a 


The matrix values for this second order system are tabulated in Table 
A-1 for three values of the constant aT, These matrix values may be 
used for any second order system with a hold which possesses the 

same aT value, making the proper allowance for the 1/, factor ap- 
pearing in equation a=?6, For example, a second Ec system with a 
pole at 10, a sampling period of „1С second, and a gain of 100 will 
possess the matrix given in Table A=1 for the aT =1 system with a 


gain of one”, 
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À-6, The System Matrix for a Type II Second Order System With A 
Hold Network, 
The transfer function of a tyre two system of interest may be 


written así 


=S T | 
GG (s) = Egg ors (a-28) 
S 


Then using the inverse Laplace transform one obtains; 
G(t) = K Га + 24) = ГЕ + 34-2)? | a(ten)} (а-29) 


Then, at the sampling instants, 


g, = K t + 2 С И (в-5)т) | n> 1 (2-30) 


The values for the matrix cf this system are given in Table 4-1. 
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TAPIE АУ} 


ELEMENT VALUES FOR THE SYSTEM TRANSFER MATRIX FOR VARIOUS SYSTEM 


la Type I Second Order System With A Hold Network 





G, G(s) 
B ex = e 5 T=1second , aT=l 
s” (s+1) 
n En = 8 (n+A ) 
0 0 50 „10653 
Y . 36788 Шу „61650 
D .76746 Aco „25895 
3 91245 259 94812 
4 96853 4.9 98091 
5 .988Д2 Des 099298 
6 4 p 9A 
fi 99843 (27; 229905 
е e TITAS 
1 ,00000 

(n+A ) cM (nea) ER 

5 „О2ёё 5 222237 
1025 e 5077 9 SPON, 
529 C189 215 09016 
25 „2135 3575 9596 
2.25 SO Kat 9852 





Ib Type I Second Order System With A Hold Network 











GG(s) -sT 
N = PES = 5 T = 2 seconds, aT = 2,0 
s“(s +1) 
B En Ms En+A 
0 0 52 226768 
1 1.13534 122 1.6819 
2 1722,09 259 1995505 
3 1.98416 СЕЈ 1592227 
A, 1.99786 4.2 15227221 
5 egg 


2.0000 





Іс Type I Second Order System with a Holá Network 


G, G(s) =ST 
B = ——— о Т = „5 SECOME „ат = ‚5 
5“ (5+1) 
2 En Rn ред 
0 О a5 0288 
ү .10653 2559 1255 
2 261237 2.5 23219 
3 (22503 Es „#727 
4 41221 9 43163 
5 244674, 389 59853 
6 245771 
7 42041 
8 48212 
9 150001) 
10 49563 


II A Type I Third Order System With A Hold Network 


G -sT 
EAS = "Luce. _ о Т = l second 
5< (5+1) (5+2) 

п En n+A En+A 
0 О A^ 05824 
1 233516 l9 2/22 
2 1.18686 250 1.47885 
3 1.67364 3.5 1,7983 
4 1.87626 
5 1,9540 
6 1.9830 
7 12599002 
oo 2 ¿C000 
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Type II Second Order System With A Hold Network 
G G(s) -sT 
RT . IS UN) , T 7 ,05 seconds 
K = 
> 
n Е 
n 
O - 0 
1 22305 
2 223125 
3 1 
4 250625 
5 228215 
6 255125 
Type II Second Order System With A Hold Network 
G(s) )( “st, 
so = dst » Т = ‚10 seconds 
K 3 
S 
n En пед 2 
0 0 #2 1925 
n 2428 ISo 3,000 
2 Boro 229 4.500 
3 2229 5 6.000 
4 6.75 
5 5.25 
n E aT + 1,50 


$ 
IL. 
N 


















01Х01 ЖЫ 
ко rT e dune 































































m pu ЖЕЕ ЕНЕ Н IEEE 





ЕЕН. 


































































PASO 
EAN A AA LELE 

CREA A LSSI A eee 
4 











Е 
ea e E e ISR 
Шиш мы ---- Шағыл 1 poseer NN 
Пана EEE 


аага ере 












аы аа а 
аа AA APA 
ES 
LLL LLL LLL EEE EE EEE EEE EEE EEE EF ! 
К ЕЕ EEE THES 






























































The Impulse Response For Three Systems. 
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AFFENDIZ В 


А МТНОр ОЕ INVERTING A LCOWER TRIANGULAR MATRIX 


The inversion of a lower triangular matrix is a necessary step 
in the evaluation of the response cf? a system by means of time do- 


main matrices, It was shown in section 2-6 that the output time 


response was given by equatíon 2-354; 


c | -|[т11-[т]+ ЕЗІ) is] ! (b-1) 
Therefore, it is usually necessary t^ calculate the inverse of 


|т]+ [e] =[a]. ҥ [А] зз written: 





| 0 
[ А] “| а," ал, O ove 





Then a matrix equation may be written)” 


А5 у = 8 = 


Thowsolution for Y is. 


|- (вінң (b-3) 
where |B| = | [4| [ 5j 77 | ч 
By the rules of matrix inversion; 


[в] - [S [4] ^. (b4) 


= [a] <o (ь-5) 


The equation b-2 may written for clarity a 


Y 
о о 


nt 











Ж 
i 
2 


Cay. Tyo ш (b-6) 
4 1 
N аа» а 
---- + = + y = ү 
a4 Уз I Y a Уз 3 
etc. 


a4 4 О 0 о 
0 222 О 0 е од 
0 233 0 
[a] = "ne 
мд о | 
5 | 
A ы а 
nn | 
a 
and 
(1/244 ) 0 C Ou EO 
О PITS O o 
ا‎ 0 0 (1/244) O oo 
em | | 
|o 0 m) 


Hence, [A] - сап be obtained by multiplication of the successive 
rows of [B] by (l/a11), (1/224), etc. Then it can be seen that 


the elements in the first cclumi of [B] are the values of Y| when 


1 
g] =| ° - )1,0,0,...! 
2 | 


The elements in the second column cf | 8] are the values of y] when 


as 


H] = {o, 1, О, О, „= 2 ani so оп, 


© 


> 
: 
) 
? 
مک ےی 
A‏ 





When Н | = 11, TR . one obtains the first column as, 


DP GL 
ЖЕСЕ 
EN! a 
а а з; 232 g 


. AU a mom Е зае 
Aa аі” m VUE E M 


This solution suggests a method of computation which allows these 


calculations to be easily carried cut, The elements of (ңі by 


(b-7) 





Уд 


which the rows of [В| аге to be multiviied, are recorded on the ex- 
treme right, This useful scheme, in the form of a table, is shown 


below for a Ax4 matrix, 





| 1 
Ж. E. ае — Ty 0 Q = 
i 0 ¿q 
a 
EN e = EE I. l 
24 1 422 
а Aas 
q —) x E p Das جوا‎ 1 0 A 
a а 
А 22 633 
а а а,. 
(el) - -42 ds Y оф bii b 5,3 1 y 
| 244 422 233 42 | а,, 


The rules for setting up this table are, 

i) To form the left-hand array enter blanks in, and to the 
right of the principal Jiagonal, Derive the remaining 
elements from the A elements as shown, 

ii) Commence the right hand array by entering units in the 


principal diagonal and zeros to the right of the diagonal 


4146 





iii) Calculate the remaining elements of [B] in succession 
by the following method. To obtain 0:3 postmultiply 


the row of the left-hand array level with b by the 


1) 
part column of [В| standing above EE Blank elements 
of the left-hand array are to be disregarded, 
iv) To find [A| multiply the rows of [B] by the scalar 


factors on the right, 


For example, the first colum cf із! is completed below; 
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As an illustration of this method cossider the inversion of the ma- 


tmix [А] given as; 


1 0 0 0 
[a] 8 - 0 0 (b-10) 
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13 are calculated as follows; 
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For a physical system it was zhcown in appendix A that invariably 
the element gg of the system matrix has a value of zero, For these 
physical systems one is usually interested in evaluating the inverse 
of [1] + [о | and therefore the matrix و‎ is simply equal to[I]. 

As an illustration consider a type I system with a hold that has 


the following transfer function 


=S T 
G, G(s) = Ое) > where Ki end T=1 
f 2 
5“ (5+1) 


It was shown in appendix A that the system matrix is. 


-— 


0 0 О C C 
23570 0 О 0 
2/075 „ed / 0 0 


[G]= | .9145 .7675 „379 © (5-14) 
.9685 .9145 ¿7673 ¿A 


367 
2554 „0085 9145 XD 
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un 


2579 


Then, one obtains for [І | + [б] » the principal diagonal containing 


element values of unity, Completing the table for inversion, one has; 


-.%79 = -= -- -= 
-.7675 -.3679  -- == == 
(9105 —.7675 -.9699 -- = 


ЕБ. 01,5--.7575 AE 
-.9884 -,9685 -,0145 -,7675 -.%79 


1 0 0 0 0 07171 
то 71 0 о 0 ОЕ 
-,6321 -.3679 1 0 0 0 1. 
-.400 -,632 -.%9 i 0 OE 
%,001 -,400 =.6321 =,3679 1 o | 1 
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Furthermore, for non-time varying systems the [G] matrix always 
contains equal values along any меге Left to lower right diagonal, 


жі 3 = 
¡Gj nesi be accomplished for only one 


Therefore, inversion of |1| + 
column as is obvious Ғгоп the example, This class of matrix called 
a diagonally invariant matrix (D,I.4.) in chapter three, lenis it- 


self to the use of a recursion formla for tre elements of the in- 


verse matrix. For the first column af the matrix H one may write; 


= = d + : УЛ ++ ooo + 5 
за а *k-1.1 "wy 221 Pio? (5-15) 
For example, in the previous solution, 
ba 7 -(231b33 + 230D32) (b-16) 
Since, the elements are invariant cùn the diagonal, one has; 
b33 =, =1 
(ь-17) 


апа b32 = Das 


Therefore, 


ba = -layı + agg Pag) = = („7675 + ,2679(-.%679) ) = -.6321 


In order to obtain accurate results, ıt is necessary to use 
either a desk calculating machine or an automatic digital computer, 
This is necessary in order to avoid the roundoff errors of multiplica- 
tion, For the usual requirements of one percent accuracy, the desk 


caleulator is completely adeguate, 























